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BROAD-SPECTRUM ELECTROMAGNETIC BACKSCATTER 


ABSTRACT 

An experimental investigation of broad -spectrum electromagnetic 
backscatter was performed. Continuous spectral responses of natural and 
manmade targets and surfaces were measured to examine the fine grain 
and gross variations with frequency of the backseattered return. Averaging 
in frequency, or effective panchromatic illumination, was investigated as 
a means of reducing the variance of the fine grain or fading variations. 

In addition, the use of multiple panchromatic bands, or polypanchromatic 
illumination, was investigated as a recognition means for the underlying 
spectral response remaining after fading effects were removed. 

The frequency variation of a discrete scatterer model was examined 
for a variety of illuminating conditions ranging from monochromatic through 
bandwidths far in excess of resolution requirements. The theoretical 
reduction in return variance was shown to be a function of the resolution 
bandwidth and illumination bandwidth ratio. 

Experimental measurements of return variance for several effective 
bandwidths were compared with the predictions of the theoretical model. 
The general behavior of decreasing variance with increasing illumination 
bandwidth agreed reasonably well; however, in all cases the experimental 
figures were larger than predicted. This effect was attributed to resonance 
variations underlying the fading distribution. The use of an approximate 
regression analysis for calculation of the experimental fading variance 
produced substatially better agreement. 

The spectral response variations remaining after the removal of 
fading affects were normalized and compared with similar response 
characteristics from the visual and infrared regions. Resonance effects 
on at least two frequency scales were noted and were observed to produce 
variations comparable to those measured in the visual and infrared regions . 
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Monochromatic and panchromatic images were produced to illus- 
trate the improvement obtainable with frequency averaging. The decreased 
scintillation of the panchromatic imagery was observed to materially aid 
the recognition of target tone and shape. Multiple panchromatic images 
were combined into a color display to illustrate the potential of poly- 
panchromatic illumination as a means of identifying spectral response 
characteristics . 

The results of this experiment graphically illustrate the improve- 
ment in return variance possible with panchromatic illumination. In 
addition, the spectral response variations remaining after the removal 
of fading effects were shown to be comparable to the variations encoun- 
tered in the visual and infrared regions. 
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CHAPTER 1 


INTRODUCTION 

An experimental investigation was made of broad- spectrum electro- 
magnetic backscatter. Continuous spectral signatures of natural and man- 
made objects and surfaces were obtained to examine the fine grain and 
gross variation with frequency of the backscattered return. Averaging in 
frequency was used to reduce the variance of the fine grain or fading 
signal and to produce panchromatic imagery. The use of multiple pan- 
chromatic bands separated in frequency was investigated as a means of 
identifying spectral signatures. In addition , multiple panchromatic 
images were combined into a color or polypanchromatic display. 

The statistical nature of radar return from complex targets or 
surfaces has long been recognized (Swerling, 1960). Only recently, 
however, has any attempt been made to utilize frequency diversity as a 
means of decreasing the variance of the return signal distribution. Al- 
though so-called "frequency agile" radar systems were developed some 
years ago, the frequency diversity of these systems was added principally 
as an anti-jamming feature. It was soon recognized that the frequency 
shift served to decorrelate the pulse-to-pulse returns and, in many cases, 
improve both detection and tracking capabilities (Birkemeier and Wallace, 
1963; Kosowsky et al. , 1963; Gustafson and As, 1964; Lind, 1968). 

Most analyses still considered the measure of improvement to be a deter- 
mination of the frequency shift required to effectively obtain independent 
pulse-to-pulse samples for various targets (Ray, 1966; Beasley and Ward, 
1968). 

Moore, Waite, and Rouse (1969) first suggested the use of 
continuous broad- spectrum illumination with a bandwidth substantially 
in excess of that required for system resolution. Operation in this 
fashion gives signal decorrelation within the resolution cell thus enabling 
one to obtain additional samples of the distribution. Moore and Waite 
(1969) in a study of Doppler scatterometer systems showed the difference 


in variance reduction between predetection and postdetection frequency 
integration. The only experimental results available for verifying the 
improvement in return distribution variance with frequency averaging are 
'.hose of Rouse (1968), and these show only the qualitative improvement 
of image quality for an acoustic modeling system. 

In Chapter 2 the improvement in return variance with frequency 
averaging predicted by the referenced methods is examined. The spectral 
.esponse data with variable bandwidths presented in Chapter 5, and the 
monochromatic’ and panchromatic imagery shown in Chapter 6 serve to 
illustrate and verify the predicted improvement with frequency averaging. 

Most of the work with "frequency agile" systems considers the 
limit of improvement to be reached when successive pulses are independent 
samples of the return distribution. Continuous broad-spectrum or pan- 
chromatic illumination can provide additional independent samples of 
the distribution within each pulse. This is particularly important for 
systems already operating at the maximum pulse sample rate (PRF) con- 
sistent with range ambiguity, such as side-looking imaging radars, 
imaging systems must provide not merely signal detection, but a reason- 
able measure of the differential scattering cross-section of each resolv- 
able element. If averaging is insufficient to reduce the fading or 
scintillation of the return, the image has a speckled or grainy 
appearance. In Chapter 3 a unified theory of broad-spectrum backscatter 
from complex targets is developed with relates resolution, system band- 
width and return variance or measurement reliability. 

The development of side-looking imaging radar has aroused 
interest in the potential of radar as a tool for geoscience applications 
(Pierson el al . , 1965; Moore, 1966; Rouse etal., 1969). This interest 
was part of the overall interest in remote sensing capabilities stimulated 
by NASA and by the potential offered by orbital platforms. Much of the 
conjecture concerning the geoscience applications of remote sensing 
assumes that a unique spectral signature can be defined for any target 
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or surface of interest (Parker and Wolff, 1965; Lyon and Vickers, 1966; 
Pardoe, 1969). This assumption has not been verified even for much more 
limited objectives. 

The investigation of spectral signatures is somewhat easier in 
shorter wavelength regions of the spectrum than microwave , for here the 
measurement and use of spectral response information has been a standard 
technique for some time. The users of the microwave or radar region of 
the spectrum have long been content with definition of returns from a few 
isolated frequencies usually identified merely with band designations. 
While a number of investigators have attempted theoretical descriptions 
of the frequency dependence of radar backscatter, virtually none of these 
has received adequate expe r cental verification {Spetner and Katz, 1960; 
Fung and Leovaris, 1968; Rouse, 1968). 

Polychromatic systems for measurements of scattering at multiple 
discrete frequencies have been used for some time (Grant and Yaplee, 
1957; Wiltse et al. , 1957; Ament et al. , 1959; Ohio State University, 
1963). However, the fundamental purpose of these systems has been 
to provide appropriate design information for particular monochromatic 
radar systems. With such a limited amount of spectral information, the 
question of signature content within the microwave region may hardly be 
approached. It is apparent that continuous spectral signatures over the 
microwave region are required both for verifying existing theory and for 
investigating possible resonance phenomena contributing to signature 
information. 

The experimental portion of this study provides the first known 
continuous spectral response measurements over a one-octave frequency 
range. In addition, the use of a color display for possible signature 
identification is illustrated by the polypanchromatic imagery presented 
in Chapter 6 . 

The design and construction of a system providing controlled 
broad- spectrum illumination was a substantial portion of the total effort. 
The system is capable of providing spectral response curves over a one- 
octave frequency range and with a variety of sampling bandwidths and 
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incidence angles. In addition, capability for the production of mono- 
chromatic and variable bandwidth panchromatic imagery is incorporated. 

A thorough knowledge of the system, calibration method and data pro- 
cessing is required for interpretation of the experimental results. A 
complete description of the system is provided in Chapter 4. 

The measurement program is divided into two sections: 1) scat- 
tering experiment, and 2) imaging experiment. The scattering experiment 
consisted of the measurement of spectral response curves for natural and 
cultural targets and surfaces. These measurements were performed over 
a range of effective bandwidths and across a frequency range of 4 GHz to 
8 GHz. The reduction in return distribution variance with increased 
effective bandwidth was compared with the theory developed in Chapter 3, 
The coarse structure (structure remaining after fading and scintillation 
effects were averaged out) of the spectral response curves across the full 
range was examined for comparative signature content. The scattering 
experiment procedure and results are presented in Chapter 5. 

The imaging experiment consisted of recording B-scan imagery 
of area extensive targets , including those targets measured in more 
detail in the scattering experiment. Monochromatic images showed the 
complex lobing or scintillation normally expected. Panchromatic images 
with sufficient bandwidth in excess of resolution requirements illustrated 
the predicted improvement. Multiple panchromatic images from widely 
different portions of the operating range were combined in a color display 
to illustrate the polypanchromatic method of spectral signature recognition. 
The imaging experiment procedure and results are presented in Chapter 6. 
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CHAPTER 2 


BACKGROUND 


The study of broad -spectrum radar is basically the study of the 
frequency dependence of electromagnetic backscatter. Studies of this 
dependence have proceeded along two widely separated lines of inves- 
tigation that might well be termed macroscopic and microscopic approaches. 
Macroscopic studies concerned with the variation of average backscatter 
power over a wide frequency range have been performed by both radar 
design engineers and scattering theorists. Most of these studies have 
been conducted to provide specific design information for monochromatic 
radar systems (Goldstein, 1946; Kerr, 1951; Wiltseetal., 1957; 

Grant and Yaplee, 1957; Ament etal., 1959; Edison etal., 1960; 

Cosgriff et al. , 1960). 

The advent of lunar radar stirred considerable interest in the 
frequency dependence of the backscattered return as a means of estimating 
the surface statistics from multifrequency data . Most of the lunar studies 
have used a Kirchhoff or physical optics theory to obtain an expression 
for the mean return involving the correlation distance and standard devia- 
tion of the surface. Curve fitting the resultant equations to the back- 
scattered data at different frequencies and incidence angles yields 
differing values for the surface statistics (Evans, 1957; Hayre, 1964). 
While this may indicate that frequency does have a "size filtering" 
effect, little can be said regarding the uniqueness of solutions obtained 
in this fashion (Fung, 1965). 

More recently the use of radar as a sensor for geoscience 
investigations has received considerable attention (Pierson et al., 1965). 
Such applications require that the radar return be correlated with the 
characteristics of the illuminated surface. The instrument parameters 
that may be varied in such an investigation are the incidence angle, 
frequency and polarization of the radar sensor (Rouse et al. , 1969). 
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Surface identification, therefore, requires a knowledge of the return 
amplitude dependence upon each of these parameters. 

The theoretical description of radar ground return depends upon 
the mathematical model used to describe the surface (Moore, 1969), The 
description of such a relatively simple surface as the sea is extremely 
difficult, even though it has relatively modest slopes and may be con- 
sidered homogeneous (Chia, 1968). To adequately describe land sur- 
faces with their greatly increased variability appears almost hopeless. 

If appears, then, that determination of radar return dependence for 
natural surfaces must depend upon empirical measurements. Theoretical 
description of the return will serve principally as an aid in interpretation 
and extrapolation of the measurements. 

The study of fine grain variations with frequency (or the micro- 
scopic approach) has been concentrated primarily in the area of tracking 
end search radar applications (Birkemeier and Wallace, 1963; Ray, 1966; 
Hansen, 1968). The received signal from any complex target or surface 
is subject to unwanted variations due to the structure of the target. The 
most common model used to describe such complex targets is one com- 
posed of a large number of independent point scatterers randomly distri- 
buted in amplitude and phase. The radar return is the phasor summation 
of contributions from each scattering center within the resolvable element. 
The variation of the received signal is normally termed "fading", although 
in some cases this term is reserved for variations in signal amplitude 
alone. The term "glint" denotes the wander of the apparent center of 
reflectivity as measured by a tracking radar. 

The random scintillation or fading of the radar return described 
above leads to errors in detection or tracking performance dependent upon 
the probability distribution of the return. It has long been recognized 
that improved performance may be realized by the integration or averaging 
of several return pulses, each a sample of the return distribution (Marcum, 
1960; Swerling, 1960). As might be expected, the improvement in per- 
formance is related to the estimation of the distribution mean and thus to 
the effective number of independent samples averaged. If the pulse -to-pulse 
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samples of the distribution are correlated, the effective number of 
independent samples may be significantly less than the number of 
pulses. Several authors have shown the improvement that may be 
obtained through the use of frequency diversity on the pulse-to-pulse 
basis to deccrrelate successive pulses and increase the effective 
number of independent samples (Kosowsky et al , , 1963; Gustafson 
and As, 1964; Beasley and Ward, 1968; Lind, 1968). 

The estimation of the distribution mean, even with the use of 
frequency agility, is still limited by the sample rate or pulse repetition 
frequency (PRF) of the system. The PRF of most systems is not deter- 
mined by these considerations hut by the unambiguous range requirement 
(Skolnik, 1962; Barton, 1964). The possible use of frequency diversity 
or broad-spectrum illumination to provide independent samples of the 
target distribution within a single pulse period has been suggested 
(Moore et al,, 1969). This technique in essence trades bandwidth for 
a better estimate of the distribution mean. A more detailed analysis of 
the relationships between bandwidth, resolution and variance of the 
sample mean is presented in Chapter 3. 

The use of the imaging radar for geoscience applications offers 
an interesting combination of user requirements. The coverage or swath 
width usually should be relatively broad, but this limits the 
pulse sample rate due to the unambiguous range requirement. At the 
same time, the geoscientist desires resolution compatible with the 
smallest features to be discriminated and a high confidence level that 
the return amplitude represents the differential scattering cross-section 
for each resolution cell. These requirements aid both detection and 
delineation of homogeneous areas. Since coverage and delineation of 
homogeneous areas are at least as important as resolution in this 
application, the use of panchromatic or broad-spectrum illumination 
offers an additional degree of freedom to the system designer. The 
ultimate aim of the geoscientist is not merely delineation of differences, 
but actual recognition. To achieve this end the addition of multiple 
panchromatic bands or poly panchroma tic illumination to determine 
spectral signature information may well be desirable. 
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In the remainder of this chapter we shall examine the theoretical 
and experimental background of panchromatic and poly panchromatic 
illumination. While the emphasis of this study is on the application of 
these techniques to imaging radar systems, virtually all of the back- 
ground material is drawn from studies conducted with other goals in 
mind. The chapter sections are further divided into panchromatic 
analyses wherein the broad-spectrum illumination is used to reduce the 
scintillation effects, and polypanchromatic analyses to investigate 
the possible information content of spectral signatures. 


2 . 1 Theoretical Analyses 

2.1.1 Panchromatic Illumination 

Numerous articles dealing with the performance of frequency 
agile radars are applicable here. The basic concept in each remains 
that of achieving pulse-to-pulse decorrelation through the use of 
frequency diversity. A sample of these papers is examined along with 
the panchromatic analysis of Moore, Waite and Rouse (1969). 

2 . 1 . 1 . 1 Birkemeier and Wallace Model 

Birkemeier and Wallace (1963) examine the improvement in 
tracking performance of monopulse and conical scan systems afforded by 
pulse-to-pulse frequency modulation. The frequency shift is obtained 
by a Linear frequency modulation with a period much longer than the pulse 
repetition period. The target model is one investigated earlier by Delano 
(1953) and Muchmore (1960) and is shown in Figure 2.1. 

The elements of the target are considered to be uniformly dis- 
tributed along the strip, each with an independent random phase angle 
uniformly distributed from zero to 2 tt , The critical frequency shift, 4fc, 
necessary to consider samples of the return signal essentially uncorrelated 
is calculated for the above model. 
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Figure 2,1 Line target model (from Biricemeier and Wallace, 1963). 

The voltage return from the target as a function of frequency is 
given by 

V(f)~ Z**cos (*¥**-+ *>) 

h 

where a = amplitude of the return from the reflecting element 
r n = distance in the range dimension from the n n reflecting 
element to a plane perpendicular to a line connecting 
the radar with the mean center of reflectivity of the 
target 

0 = independent random phase angle uniformly distributed 

n 

from zero to 2 w , 
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The autocorrelation function of the fading as a function of frequency is 
then 


R v (&f) = E[H a,cos( 


4 7T f 


ft) n. <*■» c °* ( 


f 


ft) 


( 2 . 2 ) 


where f ' = f + A f 

Since the V^'s are independent 

E [cos <P* cos p*. ] = J <C , £ f P* ' r " 1 = 7 <C 

o/i c/ £[cos<P*Stnp^] - 0 

where X n = 1 if m = n, n = 0 if m / n. 
m m 

Using the above, equation (2,2) reduces to 


R v (4 f) = f [ i L c f /J ] 

r f the frequency sweep is symmetrical such that the probability is 1/2 
.hat af is positive and 1/2 that 4f is negative, then 
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The summation of equation (2.3) may be approximated by an integral over 

the target length L noting that r * x sin © , and that the integral is an 

n n o 

even function of x, 


R.M. ^ " 

o 

ffg.J f c* a-f) 
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Let E(f) be the voltage as a function of frequency after square law 
envelope detection of the signal. Since the spectrum of E(f) is given ir;v 
the convolution of the spectrum of V(f) on itself , the autocorrelation 
coefficient of E(f) is given by the square of the autocorrelation coefficient 
of V(f). This yields 
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The decision as to how much decorrelation is required to produce 
essentially independent samples is somewhat arbitrary, however, Birkemeir 
and Wallace define the critical frequency as the first zero of the auto- 
correlation coefficient. Thus, 
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As might be expected the critical frequency change is that 
required to change by 180 degrees the phase at the target extreme and 
the phase at the reference point at the center of the target. 


2 . 1 . 1 . 2 . Ray Model 

Ray (1966) uses a slightly different model to calculate the 
frequency change required to produce essentially uncorrelated samples. 
Assume a target with a maximum dimension D', measured 
along a radius from the radar. The maximum dimension may be deter- 
mined by either the physical size of the target itself or the range 
resolution of the system. Let x be the radial distance from the nearest 
portion of the target complex to an arbitrary element of the complex, 
f the frequency transmitted by the system, and 4f a small change in f , 
For a CW transmitted signal, the phase of the echo from any 
element referred to the phase of the element nearest the radar is 



4 7T * f 
C 


(2.7) 


where c is the velocity of light. The derivative of ^ with respect to f is 


J 9 _ 4 7 TA 

Jf C 


( 2 . 8 ) 


and forms a spectrum for rate of change of phase with respect to trans- 
mitted frequency. Since the target complex is assumed to be randomly 
located scatterers with a range of x = 0 to x = D' , the result is a power 
spectrum of width AVD' /c . 

In analogy with the autocorrelation function for noise with a 
bandwidth of W given by 


&(*) 


j//. ^ tr IVt- 
Z TT W V 


(2.9) 
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Ray expresses the autocorrelation function with respect to frequency as 


e(A-f) « 


S [ h 4 7T O '_A f/c 
4 7 T O'A-f/c 


( 2 . 10 ) 


Ray assumes sufficient decorrelation to be achieved when the auto- 
correlation function falls to 0.5 which yields a required frequency 
change of 


At = 


O, (_£ C 

D' 
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( 2 . 11 ) 


If the decorrelation required is extended to the first zero of the auto- 
correlation function as done by Birkemeier and Wallace, the result is 
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c 

4 O' 
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( 2 . 12 ) 


This value is seen to be exactly one-half that obtained in the analysis 
by Birkemeier and Wallace . This difference is not due to the fact that 
Birkemeier and Wallace use the energy autocorrelation function while 
Ray uses the signal or voltage autocorrelation function since both should 
yield the same result when assuming the first zero of the autocorrelation 
function constitutes independence . The difference is due to an error by 
Ray in equation (2. 10) in which the image or negative frequency portion 
of the power spectral density is neglected. For a real spectrum of width, 
4?rD , /c, the autocorrelation function with respect to frequency should be 

dfflfj = — — n - P ‘ A (2. 

2 7T O * A f/ C 
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This yields for assumed independence at 0.5 a value of 


- 0.3 c 9 0 kAU 

Af = ■■■■ = — r Mtii (2.H) 

O D 

]f the decorrelation required is extended to the first zero of the auto- 
correlation function, Ray's result now agrees with that of Birkemeier 
and Wallace, 
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2. 1.1. 3 Moore, Waite and Rouse Model 

Moore et al. (1969) propose a system with panchromatic 
illumination where the bandwidth of the illuminating signal is in 
excess of that required for system range resolution. The excess band- 
width (that above the resolution requirement) can be used to frequency 
average the return and reduce amplitude (fading) and angle-of-arrival 
(glint) fluctuations due to phase interference effects , The effect of 
this averaging is identical to that produced by averaging independent 
samples of the return obtained by pulse-to-pulse decorrelation with 
frequency shift. The significant point is that this technique produces 
averaging within a single pulse period and hence may be used where the 
PRF is dictated by other constraints. 

Scatter from a rough surface like the ground may be expressed 
by the integral (Beckmann and Spizzichino, 1963) 
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where C - the backscattered field 

k - — the wave number expressed as a vector from dS toward 

the radar 

R - the distance from the radar to a reference point on the target 

r - a vector from the reference point to the area element dS 
f - a function of polarization, incident angle, wave number, 
and local reflection coefficient. 

Although the general character of the field is set partly by f, the detailed 
structure is almost entirely due to wave interference effects expressed by 
the phase factor. 


k ■ t = k h cos- e 


With broad-band illumination, the field observed at any point in space is 
frequency-averaged so that the fine structure of equation (2.16) is not 
present. 
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The bandwidth over which a radar must transmit to achieve the 
kind of averaging indicated by equation (2.17) depends on the fine struc 
lure in the pattern of frequency variation of the field of equation (2.13). 
The concept of required bandwidth may be illustrated by calculating it 
for the special case of a slightly rough one-dimensional surface for 
which shows (Rice, 1951) 

= ak'J ef* 2(*> (2 

' OO 
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where q = 2k sin© (0 is the incident angle) 

A = a function of distance and radar parameters 
Z = the height of the rough- surface point at location x. 

The field E is seen to be the Fourier transform in q of the surface height 
function. A power pattern may be defined as 

2 

w(t) = = w (2.i9) 

A 2 K 

where means "Fourier transform of". The power pattern is therefore the 
power spectrum of z(x). The power spectrum of a function is the Fourier 
cosine transform of its autocorrelation function, so 

W (<]) ~ j c [f* (*■)] (2.20) 


where ^ is the Fourier cosine transform symbol. Thus, the power pat- 
c 

teming is related to the surface-height autocorrelation function by a 
simple transform. 

The spacing in q (and therefore in frequency) of independent 
samples can be determined from the autocorrelation function of W, which 
is 
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As an example, consider an autocorrelation function exp f-x/L], which 
has been shown to describe surface configuration at contour map scale 
in several areas (Hayre and Moore, 1961), The expression for C{dq) is 

-(Lai/i)* 

c(a<r) = e (2.22) 


If we consider independence essentially achieved for C(Aq) = e " 2 , the 
resulting expression for frequency spacing of independent samples is 

4f= — = !_3£_ MHi (2- 23) 

7T L Slh 0 L $i* O 

which is in close agreement with the value predicted by both Birkemeier 
and Wallace and Ray, 

2.1,2 Polypanchromatic Illumination 

Polypanchroma tic is used here only to designate that the frequency 
dependence of interest is that exclusive of phase interference effects, 

2.1 .2.1 Kirchhoff Model 

The method of physical optics, or the Kirchhoff method, assumes 
that the field incident on an irregular surface is locally reflected with the 
Fresnel reflection coefficient. Stated in another way, the reflection from 
each point on the surface is considered to be identical to that which would be 
reflected by an infinite plane tangent to the surface at the point of 
observation. While this assumption does not restrict one to any particular 
range of roughness relative to the incident wavelength, it does restrict 
one to "locally flat" surfaces where the radius of curvature of the sur- 
face is large with respect to the incident wavelength. The most commonly 
used criterion for applicability is due to Brekhovskikh (1952) and is 


given by 


4 7T cof O >> A 


where f is the radius of curvature of the surface at the point of incidence 
md 9 is the local angle of incidence between the incident energy vector 
.ind the unit normal to the surface at the point of incidence. 

For a plane wave incident upon a random surface, the more or 
less standard Kirchhoff method predicts an average far-zone backscatter 
power of the form (Hagfors, 1966), 

p(e) = f JoCeksihOT*) Jr / 2 

(a cos e) z J 

where k = 2tt/A the wave number 
9 - the incident angle 
T = the distance along the surface 
J 0 = the zero order Bessel function 

X = the characteristic function associated with the surface height 
distribution 

and the integration is over the illuminated area. A proportionality con- 
stant unimportant to our discussion has been omitted from the above 
expression . 

Most theories assume a gaussian surface height distribution 
which gives 


p(0) = 


l_ 

(A cos 9) 2 
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where K = 16 ^(<T/a)^cos ^0 

0" = the standard deviation of the surface 

!’('?) = the surface autocorrelation function 

The frequency dependence of equation (2.25) is determined by 
the argument of the zero-order Bessel function and the constant K. The 
contribution to the integral of the Bessel function decreases with 
increases in either incidence angle or frequency. The factor K incorporates 
the dependence of the integral upon the relative roughness of the surface 
with respect to wavelength through the ratio (fr/U ) . As K becomes much 
larger than one / the contribution of the exponential becomes negligible 
except near the region of zero tag (r(r) close to one). The region of r(r) 
that contributes to the value of the integral is therefore wavelength 
dependent. 

The description of the wavelength dependence of the surface may 
be seen to be intimately tied to the description of the surface through the 
autocorrelation function. The most commonly assumed autocorrelation 
functions are the gaussian and exponential, although composite functions 
with wavelength dependence have likewise been used (Davies, 1954; 

Hayre and Moore, 1961; Fung and Moore, 1966), 

The various types of autocorrelation functions used to evaluate 
the integral of equation (2,22) predict a wavelength variation of A° to A 
No attempt is made to justify any one of these choices since none of them 
may serve as more than a guide to the type variation expected from natural 
targets . 

2 . 1 . 2 . 2 Spetner and Katz Random Scatterer Model 

Spetner and Katz (I960) consider the radar return from a surface 
assumed to be composed of a large number of incoherent independent 
scatterers. The normalized radar cross-section is the product of the 
density of scatterers per unit area, and the average cross-section of a 
single scatterer. At sufficiently short wavelengths , the effective 
scatterers are distinct and may be identified with the real or physical 
scatterer. 
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As wavelength increases, adjacent scatterers tend to merge and the effec- 
tive density of scatterers decreases. The density of effective scatterers 
is given by 


ff-o A < ( c ,// a o )'* 


( 2.2 6 ; 


where p - density of effective scatterers per unit area 
s 

p Q = density of real scatterers per unit area 

X = incident wavelength 
C j = constant 


The average scatterer cross section, for wavelengths small com- 
pared to the dimensions of the scatterer, is 



47 T A z 



A 4' M 0 1 I 


where p = power reflection coefficient 

A = effective area of the scatterer 


(2.27) 


For wavelengths large compared to the scatter dimensions, Rayleigh 
type scattering is assumed and the scattering cross section is given by 

C, = ~7~ ' A loh 3 e (2-28) 


where v = volume of the scatterer 
C 2 = constant 
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The transition wavelength between the application of equation (2.24) and 
(2.25) is denoted \ Q . Thus 



A ^ Aq 

A > A <i 


(2.29) 


Combining equation (2.26) and (2.27) gives four possible combinations 
for radar backscattering cross section, <j° . These are summarized in 
Table 2.1. Figure 2.2 shows in graphic form the wavelength dependence 
of a ° in the various regions. 


2 . 1 . 2 . 3 Spetner and Katz Specular Point Model 

For this model the surface is assumed to be reflective, of irregular 
shape, but continuous and with continuous derivatives. The contributions 
to the radar return are assumed to come predominately from two sources: 
the energy which is reflected from large facets of the surface oriented 
perpendicular to the line-of-sight, and the energy which is scattered 
isotropically by the small scatterers. 

The scattering coefficient is assumed to have the form 


< 7 ° = 


(2.30) 


where P ss 'P sL = scatterer densities for small and large scatterers 

respectively 

<r ls ,c 1L = scattering cross sections for small and large scatterers 
respectively 

The dependence of p is the same as shown in equation (2.26; but 

b b 

here it is valid only for large X. The average number of facet-type scat- 
terers per unit area, P sL # is the average density of points on the surface 
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TABLE 2.1 
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J ^ X 4 

x >(C 1 / P 0 ,1/2 
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o _ 1 K 

o C 1 C 2" 2 

0 * x 4 
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WAVELENGTH 


Tigure 2.2 Wavelength dependence of normalized radar cross section 
(random scatterer model) from Spetner and Katz, 1960. 



WAVELENGTH 


Figure 2.3 Wavelength dependence of normalized radar cross section 
(specular point model) from Spetner and Katz, 1960. 
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oriented normal to the incident energy. is, therefore, dependent upon 

the slope spectrum of the surface, and in particular upon the variance of 

2 

the surface slope, S . Assuming a flat-slope spectrum with a cutoff at 

_ 1 ° 2 
\ , define a wavelength dependent effective slope variance, S , related 

2 

to the measured variance , S , by 

o 


£ 


* 



i 
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A < A 2 
A > X i 


(2.31) 


The density of horizontal facets is then 


85 



A < a 2 
X > X z 


(2.32) 


The scattering cross sections for large and small scatterers are 


<v = f 


4 T[A 


(flat plats) 


(2.33) 


and 



(Rq/ UtJ k fc qtteh) (2.34) 


The transition region between the applicability of equation (2.33) and 
(2.34) will be denoted by \ 2 as in equation (2,31). Thus 
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The final expression for the normalized radar cross-section of 
an irregular reflecting surface having a gaussian slope distribution and 
a flat slope spectrum suarply cut off at V = is given by 



4 ttA* Cj (- c o f * 0/ 2 ) 

A* X l 6 

q ( - c o"t * O ) 2 Sq A £ 


/ A < A z 

(2.36) 

, A > A* 


The results for are shown in Figure 2.3. 


2 .2 Experimental Analyses 

2,2,1 Panchromatic Illumination 

While considerable has been written regarding the improvement 
in fading and glint statistics to be expected with frequency diversity, 
remarkable few measurement data have been published. 

2 . 2 , 1 . 1 . Visible Light Analogy 

Perhaps the most striking illustration of the averaging of fading 
with bandwidth may be seen in the visible region. Figure 2.4 shows 
photographs of a target illuminated first with monochromatic laser light 
and then with normal panchromatic or broad-band incoherent light. The 
scintillation or grainy appearance of the laser illuminated photograph 
well illustrates the problem of phase interference effects, while the 
smooth tones of the panchromatic photograph show the effect possible 
with frequency averaging. 

2. 2. 1.2 Philips Teleindustri Experiments 

Gustafson and As (1964) measured the comparative performance 
of a fixed-frequency radar and a so-called jumping -frequency system. 
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DE/CRES/Umv. of Kansas/ 1970 


transparency size | 5mm 

0mm 

Speckle effect illustrated by illuminating a photographic 
transparency with: 

A) Coherent (6328 A, laser) light. 

B) Incoherent (3200° Kelvin, incandescent) light. 


Figure 2.4 Effect of illuminating bandwidth in the visible region of 
the spectrum . 
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The jumping-frequency system employed frequency shifting on a pulse-to- 
pulse basis. Measurements were conducted in both search and tracking 
systems . 

The experimental results for the search mode of operation are 
summarized in Figure 2.5. The abbreviations FF and JF refer to fixed- 
frequency and jumping-frequency operation, respectively. The theoretical 
curves shown are calculated on the assumption that the fixed -frequency 
returns are completely correlated within a single scan, but totally uncor- 
related scan-to-scan while the jumping-frequency returns are assumed 
pulse-to-pulse independent (Marcum, 1960; Swerling, 1960). 

Tracking measurements on aircraft with both modes of operation 
indicated roughly a 20% increase in range to the first drop-out of track- 
ing when jumping -frequency operation was used. Figure 2.6 shows the 
i eduction of glint obtained when using jumping-frequency radar to track 
sea targets, 

2 . 2 . 1 . 3 Norden Experiments 

Kosowsky, et al. (1963) conducted angle-of-arrival and scintil- 
lation measurements for natural terrain using fixed-frequency and pulse- 
';o-pulse frequency shifting. The redistribution of the scintillation 
energy obtained by frequency-shifting is shown in Figure 2.7. The low 
frequency spectrum is substantially reduced over that obtained without 
frequency shifting, and after filtering leads to an improved estimate of 
the true angle-of-arrival in a finite observation time. 

2 .2 . 1 .4 University of Kansas Acoustic Experiments 

Rouse (1968) conducted acoustic imaging experiments of complex 
targets both with fixed -frequency and with a linear FM sweep. Post- 
detection integration was used to obtain large effective bandwidths. 

While no measurements of variance reduction were made, the resulting 
images illustrate the improvement obtainable with panchromatic illumination; 
the images of Figure 2.8 clearly show the improvement in image quality much 
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Figure 2.5 Results of experimental determination of blip-to-scan ratio 
versus range for air targets , and comparison with theoretical 
results. R = radar range; Rq = radar range giving 50% blip- 
to-scan ratio on fixed frequency. Curves 1 give the result 
of the measurements with a Jet plane (a small fighter) , 
curves 2 with a small one-engined propeller aircraft and 
curves 3 with a medium- si/ed twin-engined propeller air- 
craft (from Gustafson and As, 1964). 
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Figure 2.6 Reduction of glint using Jumping-frequency radar to track 
sea targets, as compared with fixed-frequency radar. 
Relative diatrllnition of target positions , recorded with a 
16 nun film camera and measured at intervals of 0.S milli- 
radians, is shown together with the silhouettes of the tar- 
gets, which are reproduced In their apparent actual aize in 
milliradians (from Gustafson and As, 1964). 
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Figure 2 . 



(b) Swept frequency image of two spheres 
(f= 1.5 MHz + 10%) 


8 Effect of frequency average on image of sphere targets, 
(from Rouse/ 1968). 
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1 ke that shown in the comparison of laser and panchromatic illumination 
i i the visible region. 

2.2.2 Polypanchromatic Illumination 

No one has yet produced polypanchromatic or continuous spectral 
cat a in the microwave region. A number of investigators have made poly- 
chromatic measurements at a limited number of discrete frequencies. 

Tome measurements of the effect of polarization are available at these 
limited frequencies, however, most programs have concentrated on the 
\ariation of backscutter with incidence angle. The aim of most prior 
programs incorporating polychromatic measurements has been to produce 
design data for specific monochromatic systems. 

2. 2. 2.1 Naval Research Laboratory Experiments 

NRL has operated for several years a four-frequency radar system 
mounted in a WV-2 aircraft. This system operates at 428 MHz (P-band), 
1228 MHz (L-band), 4455 MHz (C-band), and 8910 MHz (X-band). The 
antennas look forward along the fi.ght track and the angle of incidence 
•nay be varied by antenna movement. Beamwidths of the four antennas 
are 12° (P-band), 5.5° (L-band), and 5° (C- and X-band). At each 
frequency the system is capable of transmitting or receiving either a 
vertical or horizontal polarized signal, thus the full polarization matrix 
may be recorded (Ament, et al., 1959), 

The system described above has been used to record backscatter 
characteristics from a wide variety of targets, although the primary 
usage has been the measurement of sea clutter. Due to the namow beam- 
widths the variation in angle data must be collected by making multiple 
passes with different antenna tilts. While this poses little problem over 
the ocean (which may be considered homogeneous over large stretches) 
the difficulty in repeating flight lines leads to a rather gross definition 
of terrain categories. Figure 2.9 is a sample of the NRL aircraft program 
data for the categories "New Jersey Woods" and "City of Chicago". 
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ANGLE OF INCIDENCE 

(a) City of Chicago 



ANGLE OF INCIDENCE 

(b) New Jersey woods 

Figure 2.9 Naval Research Laboratory airborne backscatter 
mearurements . 
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NRL also supported a brief three-frequency ground measurement 
program (Grant and Yaplee, 1957). The operating wavelengths were 3.2 
cm, 1 .25 cm, and 8.6 mm. Measurements were recorded for a wide 
variety of terrain types and with reasonable definition. Figure 2.10 shows 
examples of these measurements. 

2 . 2 . 2 . 2 Sandia Corporation Experiments 

Sandia Corporation conducted a measurement program using a 
two- frequency (415 MHz and 3800 MHz) pulsed airborne system (Edison 
et al., 1960). Angular variations were obtained by tracking the annular 
illuminated region provided by narrow pulse modulation. Unfortunately, 
the program objective was to provide altimeter design data; hence, the 
technique used restricted the angular range to near vertical incidence 
less than 30° from the nadir). The targets measured are reasonably well 
documented and samples of the measurements are shown in Figure 2.11. 

2 . 2 . 2 . 3 Ohio State University Experiments 

Prooably the most carefully controlled and repeatable series 
measurements of natural surfaces has been conducted by the Ohio State 
Jni versify Antenna Laboratory (Cosgriff et al , , 1960). A three -frequency 
truck-mounted system operating at 10 GHz (X-band), 15,5 GHz (Ku-band) , 
and 35 GHz (Ka-band) was used. The antennas were mounted on an 
elevated boom and were adjustable over a range of incident angles. The 
measurements were taken by moving the truck along an assumed homo- 
geneous stretch of terrain while recording. 

The principal disadvantage with this system lies in the relatively 
small illuminated area (approximately 1 meter in diameter) due to the low 
altitude of the antennas. This restricts measurements to small scale 
terrain features such as grass, concrete roads, etc. The apparent accuracy 
and repeatability of those measurements offers a valuable insight 
into the behavior of radar backscatter from natural terrain. Figure 2.12 
shows examples of Ohio State University data. 
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Figure 2.10 Naval Research Laboratory ground-based backscatter 
measurements (after Grant and Yaplee, 1957). 
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(a) Farmland near Cameron, Missouri. 



Figure 2 .11 Sandia Corporation airborne backscatter measurements 
(after Edison et al., 1960). 
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7 . 2 . 2 , 4 Waterways Experiment Station Experiments 

The U.S. Army Waterways Experiment Station operates a four- 
frequency system mounted on a 50-foot radius wooden arch. The antennas 
v travel along tracks mounted inside the arch and may view an approximately 

< T -foot diameter target sample at the center of the arch over an angular 
range of 210 degrees. The system frequencies are 297 MHz (P-band), 

5 879 MHz (C-band), 9375 MHz (X-band), and 35, 543 MHz (Ka-band). 

Measurements conducted with this system have been principally 
concerned with determining the amount of penetration into soil and the 
effect of moisture content upon the backscattered return (Lundien, 1966). 
,r, hese measurements are made on carefully prepared soil samples. A 
limited amount of data from more natural targets has been taken and a 
sample of these results is shown in Figure 2.13. 
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Figure 2.13 Waterways Experiment Station ground-based backscatter 
measurements (after Lundien, 1966). 
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CHAPTER 3 


VARIANCE REDUCTION OF RADAR BACKSCATTER WITH 
v BROAD-SPECTRUM ILLUMINATION 

I 

l ■ 

In this chapter we shall investigate the effect of broad- spectrum 
illumination upon the fading statistics of a target model assumed to 
represent a distributed target. It is assumed that the targets of interest 
are sufficiently complex that they may be represented by an equivalent 
set of point scatterers . The behavior and, in particular, the variance 
of the signal return from the postulated model is then examined for a 
variety of transmitted signals both monochromatic and broad -spectrum . 


3 . 1 Target Model 

We shall assume that the targets of interest are sufficiently 
complex that they may be represented by a random collection of dis- 
crete independent scatterers. Each scatterer reradiates isotropically 
with an amplitude proportional to its effective radar cross-section. 

The amplitude distribution is random and is assumed to have a gaussian 
distribution with zero mean and variance, (F - . For a given target area 
or class the average density of scatterers per unit area,!", is constant, 
i.e scatterers are uniformly distributed within the target. 

Figure 3.1 shows the target model illuminated by a radar system. 
Angular resolution of the target is provided by the antenna beam width, p . 

The antenna gain is assumed constant over the illuminated area and 
zero elsewhere. Range resolution is obtained by modulation of the trans- 
mitted waveform. The dimensions of the target may be either larger or 
smaller than the corresponding dimenions of the radar system resolution cell. 
For instance, the resolution cell of an imaging radar system may occupy 
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only a small portion of a large homogeneous field. At the other extreme, 
the target for a search radar system is frequently an aircraft with an 
area or volume much less than that of the system resolution cell. 


Radar 



Figure 3.1 Discrete scatterer model of target area. 


Within the target complex each scattering element is characterized 
by its return signal amplitude and its range from the radar. In all 
instances it is assumed that the target range dimension, D, is much less 
than the minimum range of the system such that the return signal amplitude 
of the scattering elements may be considered independent of range. 

The signal returned to the radar from a single scattering element 
is 


S, <*) = Q h S(*-*») 


(3.1) 


where s(t) = transmitted waveform 

t = round-trip time delay corresponding to the scatterer 

range 

a n = scattering amplitude adjusted for system parameters. 
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The local return from the target is the summation of the returns from 
all scattering elements located within the illuminated area. Thus, 

N 

Sf(t> = z. (3.2) 

h = I 


The delayed transmitter signal may be expressed as the convolution of 
the transmitted waveform with an impulse function. 

<f r*-A) <!■*•> (3-3) 

— oo 


where 


J (Jt 


(o if X * 

| / . f X = X h 


The total received signal may now be expressed as 

N *o 

S h (A) = ]T Q* J S(*h) < 3 - 4 ) 

h ~ I ' <*> 

Interchanging the order of integration and summation, we have 

50 A/ 

S.OO = J S(*«) Z a - S (*-■**) (3.5) 

*"oO H'Zl 

The return signal is, thus, the convolution of the transmitted signal with 
what might be considered the impulse response of the target model. We 
may express equation (3.5) in a simpler fashion by 


s h a) - rw <9 cU) 


(3.6) 
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where the symbol, ® , denotes convolution and c(t) is the pulse response 
of the target model given by 

N 

cf*)= ]T Q* (3.7) 

h c I 


The frequency response, C(f), of the target model is given by 
the Fourier transform of the impulse response. Thus, 

C({) = J H f c it 

~ oO 


N 

- r 


2 7r f 

e 


* a / 


(3.8) 


C(f) may be characterized as a vector sum with random amplitude and 
phase. For any frequency f >>1/0, the phase distribution will be 
effectively uniform over the basic phase interval 0 to 2 te. By taking 
the minimum range to always be significantly greater than zero, it is 
assumed that the effective scattering amplitudes are not only independent 
of range, but that they meet the requirements of the Central Limit Theorem, 
i.e the variance of any one ^erm is negligible compared to the variance 
of the sum . 

With the above conditions the random vector sum is recognizable 
as a complex gaussian variable with a Rayleigh distributed envelope and 
a uniformly distributed phase. Thus, if 

B = I c (f) I 


we have the density functions 


/a(l) = 


2 



(3.9) 
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f. (e) = 


nr 


0505 l 7T 
e /f e 


(3.10) 


where 


- i-E f[<l ‘if C-. 


h = l 


h = I 


(3.11) 


The mean, mean square and variance of the Rayleigh distribution are given 
by 


ytj = e[i] = 

(3.12) 

£ [i l ] = 2 * * 

(3.13) 

<?* = £[(*-?>*] = 

(3.14) 


Thus, the magnitude of the spectral response of the target model is given 
by a Rayleigh distribution at any frequency. 

The frequency response autocorrelation function of the target: 
model may be obtained from the expression developed for C(f), thus 


R c (V) = £■[ c(f) ] 


= £ 


h = / 


~ I 


= £ [ E E>. ^ it”™'-*'* ] 

h sr | tnz. I 


(3.15) 
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Since the amplitudes of the individual scattering elements are statisti- 
cally independent, the expected value of the double summation will be 
zero except when n = m. Therefore, 


p™ r * ~ j 2 TT )) £ h *1 

Red) = E e/ J (3.16) 

h ^ J 

The strength and position, i.e. amplitude and phase, of the individual 
scatterirg elements are statistically independent. So 

R c (v)= II £[*•!]£[&*”'*''] < 3 - 17 > 


We further assume that not only are the scattering elements uniformly 
distributed, but that the distribution of time delay is likewise uniform. 






For the case of a homogeneous area target of greater extent than the system 
resolution, this involves approximating the angular sector of Figure 3.1 
by parallel lines over the range interval of interest. Taking the target 
center as the zero time reference, we have 

s &*"**'■ 

- mo 

5 ^ / -/ Z IT Vj* 

= / T e/ tU * 


SI h 7 r T c 
TV! c 


(3.18) 
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(3.19) 
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We now have 


R c fv) 



ft* 7 T V> T c 
7T V> 7^ 


(3.20) 


Recalling that T = 2D/c, we have in terms of target size 
c 


/?*(*> 



,f / n £ 77 ^ p/c 

2 7 rv £>/c 


(3.21) 


The spacing between independent frequency samples of the model 
may be determined from the autocorrelation function of equation (3.21). 

If we assume independence to be achieved at the first zero of the auto- 
correlation function, the required frequency spacing, Af, is given by 



tJ-Q 

D 


MH 2 


(3.22) 


The above result agrees with those shown earlier in Chapter 2, 
which is not surprising as the model is basically the same. We may now 
examine the effect of this model upon the return signal. If we transmit 
a signal, s(t), with corresponding spectrum, S(f), the received signal 
from the target model will be given by 

*2 if { Jt 

$,.<*) = s(jt) © c (*) = f s(f) c(f) y f (3.23) 

* 04 
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where c(t) and C(f) are the impulse response and corresponding spectrum 
of the target model. The total received energy from the target is given by 




= f /f,. <■*)/*<// = flS(f>l* lc(f)l* J f 


(3.24) 


The received energy, W .is itself a random variable that will change as 
we examine different sample functions of the target process. The expected 
value of the received energy is 


E[w t ] - E[flS(f)\ l lC(f)l l Jf 

~ GO 

= E[lC(f)\ t /°l Jf 

" GO 


(3.25) 


The integral portion of equation (3.25) is recognized as the transmitted 
signal energy, W , while the mean square value of C(f) was previously 
determined to be 

£ [ 1C (-f)l*] = = E ' G* 

h = » 

The expected value of the total received energy, then, is given by 

E[Wy] = (T< l W t (3.26) 


Of perhaps even greater interest is the variance of the received 
energy about its expected value. This is given by 


< = e[wS ] - £'[*>] 


(3.27) 
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2 

We may express W f as 

w* = f IS(f)l t IC(f)l t Jf fls(f')l 2 lc(f')l*Jf' 

~ OO “ « 

CO «c 

=yY ism^icmi* <>f jv 

--O (3.28) 


where f' = f - v 


The expected value of W z is 

r 

60 cm3 


E[^J - //E[lC(f)l i IC(f-»)l l ISWI l lE(f-i)l z df </v 




oo «0 

= Jet tc(f)i f -ic«-v)i*][ 


eo 


(3.29) 


The second integral is recognized as the autocorrelation function of the 
transmitted signal energy spectral density, K (v) . Thus, 


E[W h z ] ~ f E[<C(f)l z IC(f-V)l i k s (T>) dD 


OO 


(3.30) 


The target model spectrum, C(f), was shown to be a complex gaussian 
variable with zero mean and uniformly distributed phase. For this type 
variable it can be shown in general that (see Appendix A) 




(3.31) 
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Applying the above relationship to the expected value portion of equation 
(3.31) yields 


£[/crf)/ ! /C(Y-v)/*J = E*[lC(f)l*J t£ ! fcWc*, 

where E*[ lC(f)l l ] = <7 £ 4 

4 r HtULl 

and £ i [c({)C»(f-\»] = L ZttMD/c 


We now have for the mean square 

o a 

e[w* ] = f [<r t A + RUv>] KtO) ^ 

- *0 

The integral of K g (v) is equal to the square of the transmitted signal 
energy. Thus,. 

oo 

E />„*] • 4 *• <*> 

“CO 

The square of the mean is 

E e [w h ] = <r t 4 iv/ 


The variance of the return energy is then given by 

(Jw h - f°Rc Cv) H s f»J J v 

-oo 


-V)J 


2 


(3.32) 


(3.33) 


(3.34) 


(3.35) 
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The variance of the total received signal energy is / thus, a 
function of the autocorrelation functions of the target frequency spectrum 
and the signal energy spectral density. 


3.2 Monochromatic Illumination 


True monochromatic illumination implies an infinite duration 
time signal and cannot, of course, provide range resolution of the 
target. This condition might be approached in practice by a system 
sensing a self-bounded target, such as an aircraft, with a single 
frequency pulse of spatial extent greatly exceeding that of the target. 

Since the monochromatic signal has infinite energy, it is necessary to 
use the power spectral density rather than the energy spectral density. 

We shall assume a signal of the form 

soo = < 3 - 36 > 

The spectrum of this signal is an impulse and is given by 


S(f) = 


(3.37) 


The power spectral density is likewise an impulse located at the same 
frequency 


lS(f)l* = < 


(3.38) 


The autocorrelation function of the power spectral density, K g (v) is thus, 

K,C») = / lS(f)l i IS(f-v)l t c/v 


= J° S(f-f.) s(-f-f.-v) Jv 

— Op 

= SCO) 


(3.39) 
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The variance of the average power returned by the monochro- 
matic signal is given by equation (3.35) as 


2 ^ 1 ) o/c ) 


- »o 


= <rs f 


«o 


(e 7T \f d/c) 2 


- 


<f (o) <J^ 


(3.40) 


Recall from equation (3 . 34) that 

£ IVJ * 0c* Wjr 

where W is the transmitted signal energy or, in this case, the average 
s 

power transmitted 

OO 

W^f * /K(0/ ! </f = f ls(*)l‘ 

= / 




(3.41) 


Thus, for the monochromatic signal the expected value of the average 
return power is 


£[WJ = <Tc* 


(3.42) 


while 
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We note that the variance of the return is simply the square of 
he expected value. This is not surprising since the point scatterer 
model used is known to give rise to the Rayleigh distribution when mono 
chromatically illuminated. The Rayleigh power or exponential distribu- 
tion is characterized by equal mean and standard deviation, thus the 
variance of return power is equal to the square of the mean. 


3 . 3 Broad-6 pectrum Illumination 
3.31 Rectangular Spectrum 

For the case of broad-spectrum illumination we shall first con- 
sider the case of a broad-bar d rectangular spectrum of the form shown 
in Figure 3.2. 



Figure 3.2 Rectangular frequency spectrum. 

The total energy transmitted is normalized to the average power 
value obtained with monochromatic illumination. This will facilitate 
later comparison with the monochromatic result ard that of other types 
of illumination. The total energy for this signal is given by 


p*/i 

W f = / l$(f)l'Jf = J A*c/f 

/i f*/l 


= A l f* 


(3.43) 


Normalissing the above result to one requires that 



The signal energy spectral density and its autocorrelation function 
are shown in Figure 3.3. 


|S(f)|*i 
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(a) Energy spectral density 



(b) Autocorrelation function of 
energy spectral density 


Figure 3 . 3 Rectangular spectrum . 
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We have, then, forK c (v) 


H s (\» = /lS(f)l z IS(f-y»l z Jf 


- oo 


I r IV/ i . r 

--T 0-xJ 


(3.44) 


The variance of the return energy is now given by 
2 


- /*.'<« t f '* x J *'■' 


-f* 
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(7 c 4 ( Zjnf^D/c) | 

f M -/ fzn-VD/eJ* L ~ f *. J 


(3.45) 


Making the following change of variables 


= Z 7 T \) 0 /c /?= 2 ^ P/ C 


we have 


TTft 

e 2<7c 4 /• [■ _ X 1 y 

V * IFff - 7 “TT" L 


(3.46) 


Not too surprisingly, the result turns out to be dependent only 
upon the signal bandwidth-target size product. Since the time extent of 
the target, 2D/c, is inversely related to the resolution bandwidth of the 
system we might also consider R to represent the ratio of transmitted 
signal bandwidth to system resolution bandwidth. This ratio will, of 
course, depend upon precisely how the system resolution is defined. 
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While the integration indicated by equation (3.46) cannot be 
performed directly, it is relatively simple to evaluate by means of 
rational approximation (see Appendix B). The results of the approximate 
integration are shown in Figure 3.4, 

The variance of the return energy as shown in Figure 3.4 is 

4 

normalized with respect to a . This may be interpreted as a normali- 

c 

zation with respect to the variance of the monochromatic result, or the 
quantity ^ 2 /ct^ may be considered the square of the ratio of the 
standard deviation and mean of the return energy distribution. 

<. c trJ, s <?wV = \ ! 

(tf (3.47) 


It is interesting to note that the addition of signal bandwidth 
causes an appreciable reduction in return variance even for R equal to 
or less than one. For any practical system this means one should expect 
a return distribution with slightly less variance than the Rayleigh distri- 
bution which corresponds to the monochromatic example. For instance, 
the rectangular spectrum considered here will, with linear phase varia- 
tion, give a ™ 0 - ^ type time waveform. If we consider the system resolu- 

^ sin x 

tion to be defined by the separation between the first zeros of the 

X 

function the resolution bandwidth will correspond to R = 1 and the 
variance of the return energy will be 0.656 that of the Rayleigh distribution. 

3.3 2 ' ^x x Spectrum 

e j H y 

The spectrum is particularly interesting as its transform 

in the time domain is a perfectly rectangular pulse. With this pulse 
adjacent range resolution cells are completely independent as there is 
no overlapping illumination on the leading and trailing edges. We shall 
assume a spectrum of the form shown in Figure 3.5. 
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Figure 3.4 Variance reduction with rectangular spectrum 
panchromatic illumination . 
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Figure 3 . 5 


sin x 
x 


frequency spectrum . 


Again, the total energy transmitted is normalized to facilitate 
comparison with other types of illumination. 


W s = f \$(f) l* J f 

— OO 





</f 



( 3 . 48 ) 


Normalization yields a spectral amplitude of 

A * 


The autocorrelation function of the signal energy spectral 
density is given by 
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K s (v) = /lS(f)l Z IS(f~i)l Z J V 

— *0 

= /imrei”" Jr 

- GO 

(3.49) 

The time domain representation is somewhat easier to work with for this 
type signal, thus, the second integral is used. The time domain function, 
its autocorrelation function, and the square of the autocorrelation function 
are shown in Figure 3.6. 
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= -4- f Xr _ BVjXp I 

‘-('27TVJ* (2rrV)3 J Xf- ~ v - (3 . 50) 

The variance of the return energy is now given by 
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(3.51) 
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(b) Time domain autocorrelation function 



(c) Square of autocorrelation function 
Figure 3,6 — spectrum. 
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Making the following change of variables 


T c = ■ ^ ~ — Qi\ d -Y - 77" V T c 
c 


we have 



The similarity between the above result and that obtained 
previously for the rectangular spectrum is obvious. Again the parameter, 
K, is subject to various interpretations. It may as before be considered 
a signal bandwidth -target size product, a ratio of signal bandwidth and 
resolution bandwidth or a ratio of target time extent and system time 
resolution. The precise nature of the parameter, K, will again depend 
upon how the resolution bandwidth is defined. 

The integral in equation (3.53) is likewise intractable, however, 

2 x 

a series expansion of the sin(-^r) term results in a form easily integrated 
by computer. The results of the approximate integration are shown in 
Figure 3.7. 
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Figure 3.7 
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Variance reduction with spectrum panchromatic 

illumination . x 
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3 . 4 Postdetection Averaging of Independent Resolution Cells 

Figures 3.4 and 3.7 illustrate the decrease in variance of the 
tsnergy distribution obtainable with panchromatic illumination (excess 
range bandwidth). Let us compare this with the improvement that may 
be obtained with a system that processes bandwidth for maximum reso- 
lution and then subsequently averages the return from different resolu- 
i ion cells . 

We shall assume that the return distribution for a maximum 
resolution system is the same as that obtained with monochromatic 
illumination (Rayleigh distributed). As shown previously, the return 
energy distribution with monochromatic is given by 


/i (v) 


/ 

2 <7 C 2 


- 




(3.54) 


The above distribution is recognized as the Rayleigh powe,. or 
exponential distribution and also as a Chi-square distribution with 
two degrees of freedom, for 


/2(X Z ) 


(_ 

(z <r c z J (»/z -0/ 
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<Vz -o 


- x Vz ^ 1 

e 


(3.55) 


where n is the number of degrees of freedom. The mean value and 
variance of the Chi-square distribution are given by 


h) e q /i 



l /ah / h ce 


zh<r + 


The return from each resolution cell of an assumed homogeneous 
target will be characterized by the above distribution. Averaging of 
these returns involves the addition of random variables each distributed 
with a Chi-square distribution with two degrees of freedom. If we 
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assume each distribution to come from a totally disjoint segment of 
the target area the distributions will be statistically independent. The 
additive property of the Chi-square distribution says that the sum of 
independent Chi-square distributions will yield a Chi-square distribu- 
tion with the number of degrees of freedom equal to the sum of the degrees 
of freedom of the individual distributions , 

From the above we may calculate the reduction in normalized 
variance as more resolution cells are averaged, by 

g h <T C 4 _ _ 2 _ 

h z <TS ~ * (3a5 



Since each independent sample of the return distribution has two degrees 
of freedom, the number of independent samples averaged is one-half the 
number of degrees of freedom. Figure 3.8 shows the linear curve of 
normalized variance as a function of degrees of freedom and independent 
samples. 


3 . 5 Comparison of Averaging Methods 

Figure 3.9 shows the variance reduction obtainable with the 
methods investigated. Again, note that each result is normalized to 
that obtained with monochromatic illumination and that the total energies 
transmitted are equal. 

At first sight it appears that the panchromatic illumination 
methods actually produce a slightly better result than that produced by 
postdetection averaging of independent resolution cells. However, this 
appearance is due to a number of simplifying assumptions made during 
the course of the development. 

First, consider the difference between the rectangular spectrum 
and postdetectlon averaging curves. As pointed out earlier with regard 
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Figure 3.9 Comparison of averaging methods. 
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to the rectangular spectrum, the addition of sufficient bandwidth to 
merely resolve a resolution cell significantly reduces the normalized 
variance and causes a departure from the Rayleigh result obtained with 
monochromatic illumination. For instance, assuming the resolution 
bandwidth to be given by the reciprocal of the target time extent, 2D/c, 
the bandwidth -target size parameter, R, must equal one in order to 
resolve the target. This, in turn, says that the return from each reso- 
lution cell is characterized by a distribution with less variance than 
that of the Rayleigh distribution. Thus, the postdetection averaging 
estimate is actually somewhat higher than the result one would expect 
from an actual system. 

On the other hand, the rectangular spectrum is known to yield 
a - in - time signal of infinite duration. Thus, consideration of statis- 

X 

tically independent resolution cells actually involves an approximation 
as to how much overlapping illumination may be permitted while still 
assuming essentially zero correlation between cells. 

While these points are minor, they do account for the discrepan- 
cies between the two curves. In addition, one is led to the conclusion 
that the curve for the rectangular spectrum is probably the better approxi- 
mation to a physically realizable system. If this is the case, we may 
expect the result of averaging over a relatively narrow panchromatic band 
or averaging a small number of resolution cells to produce slightly less 
variance than that normally predicted by consideration of independent 
sample averages. 

Let us now consider the differences noted in the curve for the 

sin x 

— spectrum. This spectrum yields a perfectly rectangular time 

signal and, thus there is no problem in considering disjoint independent 
resolution cells. However, this spectrum has the rather undesirable 
characteristic of requiring an infinite system bandwidth. 

c|n v 

In comparing the rectangular and— 1 results we are again 

XT 

faced with the problem of defining resolution. If we define resolution 
in the time domain, this presents no problem for the S11 ^* - spectrum 
gives a rectangular pulse. The rectangular spectrum, however, produces 
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a “Y 3 ^ pulse in the time domain and some approximation must be made 

regarding what constitutes resolution. If we consider the resolution to 

be defined by the region between the first zeros of the function, 

the resolution bandwidth is given by 1/2 t f . For the same time domain 

sin x 

resolution, however, the central lobe of the spectrum will have 

twice this width or l/t r . Conversely, if we match central lobe and 
rectangle in the frequency domain we obtain a two-to-one difference in 
the time domain. 

sin x 

Applying this two-to-one difference to the result for the — — 
spectrum we have the result shown in Figure 3.10. Of course, defining 
resolution or bandwidth as extending to the first zeros of the 

A 

function is arbitrary and the difference will be correspondingly less if 
the 3 db points are selected for this definition. 

For a practical system confined to finite time and bandwidth , 
the actual results will lie somewhere between the curves shown in 
Figure 3.10. We note that each of the curves of Figure 3. 10 approach 
a straight line with the same slope. Considering the resolution band- 
width to be 



c 

2 6 


and, likewise to be the bandwidth necessary to achieve independent 
samples, we may calculate the frequency separation between independent 
samples as 


A f = f , = -5- MHa (3.57) 

This is a familiar result and merely states again the equivalence 
between what is considered an independent sample and the bandwidth 
required to resolve the sample. Again, this result is dependent upon 
the definition used for resolution. 
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From the preceding development we may conclude that the means 
of frequency averaging is immaterial so long as energy is conserved. 

With this conclusion the question arises as to what advantage is gained 
by using panchromatic illumination. The answer to this lies in the actual 
means of system implementation. While the result is essentially the 
same for a maximum resolution system with post-detection averaging and 
a panchromatic system, the point at which the averaging is performed may 
be significantly different. 

To realize any advantage from the panchromatic technique the 
averaging operation must be performed prior to the final video output or 
recording. This may yield a substantial savings in the bandwidth 
necessary in the main receiver amplification stage and for imaging systems 
in the bandwidth of the recording system . Since the recording bandwidth of 
high resolution systems is one of the prime constraints this offers the 
possibilities of a significant savings. 
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CHAPTER 4 


EXPERIMENTAL INSTRUMENTATION 


The system used to produce the measurement data for this study 
was designed and built as a portion of the research program. The design 
o: the measurement system was considerably influenced by the desire to 
produce a flexible instrument capable of a wide variety of broad -spectrum 
measurements rather than one restricted to the limited goals of this 
specific program. The choice of subassemblies for the system was 
determined as much by this desire for flexibility as by the normal con- 
straint of cost. 

To fully investigate the capabilities of broad-spectrum illumina- 
tion or "polypanchromatic" radar requires an illuminating source capable 
of extreme broad -bandwidth transmission. This is not in principle a 
difficult task since sub-nanosecond pulse-width radars have been built 
and operated by numerous investigators (Micronotes, 1967). Techniques 
tor pulse compression to increase bandwidth are likewise common in 
today's operational systems. If the spectrum of the scattered broad- 
bandwidth return is to be processed for spectral information content it is 
also necessary that the system have multiple receiving channels, each 
responsive to a different bandwidth of the illuminating spectrum. 

The above requirements for broad -bandwidth and multiple receiv- 
ing channels, while within the state-of-the-art, were well outside the 
budget for this program. It was possible, however, through the use of a 
so-called "slow-sweep" mode of operation to build a flexible system for 
investigation of a wide range of illuminating conditions. 

In line with the desire for maximum flexibility the entire system 
is mobile including its own source of power furnished by a diesel motor- 
generator set. The antenna is mounted on a 40-foot hydraulically operated 
boom and may be positioned for any incidence angle or azimuth setting by 
a remote controlled antenna positioner. The complete system is shown 
in Illustration 4.1. 
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4 . 1 System Function 


In the "slow~sweep n mode of operation the transmitter frequency 
is continuously swept across the full operational bandwidth at a relatively 
slow rate. This continuously swept signal i/s pulse modulated at a pulse 
repetition rate much faster than the sweep r/a te, thus the transmitter out- 
put is a succession of pulses each at a slightly different carrier frequency. 
By suitable choice of the relative repetition rates and the pulse width, 
the spectra of the pulses are overlapped ^giving an almost uniform output 
spectrum over the sweep range. This technique is illustrated in Figure 
4.1. / 

The average amplitude of the/return pulses is, in essence, a 
sampled data record of the target spectral response across the swept 
band. The sampling window spectrasl width is that of the individual 
amplitude modulated pulse. Postdetection averaging of successive 
pulses enables one to increase tbye effective bandwidth by small incre- 
ments up to the full sweep width/ The fine control of effective bandwidth 
permitted by such averaging in frequency permits investigation of the 
effect of bandwidth or panchromatic illumination upon the fading return 
of complex targets . 

The need for multiple receiving channels to investigate possible 
spectral signature identification is met by recording B-scan images of 
area extensive targets for several panchromatic bands. The resultant 
images are then optically or electronically combined to produce false 
color "polypanchromatic" images (Dalke, 1969). 


4.2 System Description 

The block diagram of the measurement system is shown in 
figure 4.2. The system is composed of six main sections: 1) trans- 
mitter, 2) receiver, 3) local oscillator, 4) digital controller, 5) scatter- 
ing coefficient data processor, and 6) imaging data processor. 



(a) Sweep modulation. 
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(b) Pulse modulation. 






(d) Overlapping spectra due to slow 
change of carrier frequency . 


Figure 4.1 Slow-sweep method of obtaining broad-bandwidth. 
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Figure 4.2 System block diagram 










































4.2,1 Transmitter Section 


The transmitter section consists of a general-purpose sweep 
frequency oscillator, signal splitter, PIN diode leveling loop, micro- 
wave modulation switch, power amplifier, microwave noise blocking 
switch, duplexer, and antenna. 

4 . 2 . 1 . 1 Sweep Frequency Oscillator 

This oscillator is a commercial Hewlett-Packard Model 8690A 
sweep oscillator with a Model 8693B plug-in module that covers the 
range of 4-8 GHz. The sweep output may be externally triggered to 
sweep between two variable frequency settings within the 4-8 GHz 
range. External trigger outputs are available synchronous with the start 
and stop of the linear frequency sweep as well as a blanking pulse 
coincident with the generator retrace and stabilization interval. Sweep 
time is continuously variable between 0.01 and 100 seconds, and the 
minimum power output across the full range is 15 milliwatts when leveled. 
The PIN diode leveling components originally mounted in this unit were 
removed for external use after division of the output signal. A sweep 
reference ramp voltage is available from the unit and is used to drive a 
swept frequency filter in the local oscillator section. 

In normal operation the sweep is externally triggered by the 
digital control unit to sweep from 4-8 GHz with a sweep time of 0.01 
seconds. The blanking pulse output is used in the digital control unit 
to suppress spurious trigger signals during the retrace and stabilization 
interval. 

4. 2. 1.2 Signal Splitter 

The output from the sweep oscillator is divided to both drive the 
transmitter power amplifier and to provide a local oscillator reference 
signal. This division is accomplished with a -10 db parallel line direc- 
tional coupler, Narda Model 3004-10 . The -10 db output is used for the 
transmitter drive while the straight through output is used for the local 
oscillator reference signal. 
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4 . 2 . 1 . 3 PIN Diode Leveling Loop 

The transmitter line power is leveled using the PIN modulator 
and automatic level control amplifier originally used to internally level 
the sweep oscillator. The leveling signal is obtained by sampling the 
transmitter line with a Narda maximally-flat -2 0 db coupler. Model 
3044B-20, and a Hewlett-Packard Model 423A detector , This arrange- 
ment reduces the transmitter line power variation at the sampling point 
to within + 0.4 db across the full 4-8 GHz band. 

4 .2 . 1 .4 Microwave Modulation Switch 

The swept frequency transmitter signal is pulse modulated by 
means of a Hewlett-Packard Model 3570 microwave switch. This unit 
uses two PIN diodes in a grounded cathode configuration and provides 
a 1 .5 db maximum insertion loss and a 38 db minimum isolation across 
the 4-8 GHz band. In the actual system two of these units are used in 
series with parallel switching to give sufficient isolation in the off 
condition. The switching speed of the series combination is approxi- 
mately 12 nanoseconds off-to-on and 8 nanoseconds on-to-off. With 
these speeds it is possible to operate the system with a pulse width 
of as low as 20 nanoseconds with no appreciable loss of peak power 
output. 

The relatively fast switching speeds used in the system make 
adjustment of relative switching time rather critical since the variation 
in propagation time and switch delay is significant in comparison to rise 
time and pulse width. In order to precisely adjust switching times a 
switch drive circuit was designed to permit independent delay of the 
leading and trailing edges of the driving pulse. The circuit diagram of 
the switch driver is shown in Figure 4.3. 

4. 2. 1.5 

While the peak power required to meet the system design goal 
of 2 0 watts is quite modest, the need for 4-8 GHz operation virtually 
dictates the use of a travelling wave tube (TWT) power amplifier, A 
grid-modulated TWT would seem to be called for, howevever, the need for 
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+5v 



AU resistances in ohms 

Diodes: 1N4009 Transistors: NPN 2N4420 
Cj and C 2 selected for required delay 

Input: -5 volts: switch OFF to RF 
0 volts: switch ON to RF 


Figure 4.3 Driver circuit for grounded cathode microwave PIN diode switch. 
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-5 v 


All resistances in ohms 

Diodes: 1N4009 Transistors: NPN 2N4420; PNP 2N4423 
C 3 and C 4 selected for required delay 

Input: *5 volts: switch ON to RF 

0 volts: switch OFF to RF 


Figure 4.4 Driver circuit for grounded anode microwave PIN diode switch. 
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extremely fast rise time and relatively low cost led to specification of a 
standard commercial unit with 20 watt CW capability. The CW operation 
required that external modulation be provided and further necessitated an 
extremely high duplexer isolation to keep the TWT amplifier noise from 
the receiver input during reception, The amplifier used is an Alto 
Scientific Company Model 20C 4. 0-8.0 Q60, This unit has a 20 watt CW 
rating with a 60 db small signal gain and a maximum noise figure of 35 db. 

The TWT amplifier is operated with PIN diode switches at both 
the input and output. In order to reduce the VSWR and prevent possible 
damage to the amplifier, load isolators are used at both the input and 
output of the TWT, These units are Melabs Model H-317-264-T internally 
terminated ferrite circulators with a maximum insertion loss of 0.4 db 
and a minimum isolation of 18 db across the 4-8 GHz band, 

4 . 2 . 1 . 6 Microwave Noise Blocking Switch 

The CW operation of the relatively noisy TWT power amplifier 
makes it necessary to provide an exceptionally high transmitter-to- 
receiver isolation. This problem is further complicated by the wide 
bandwidth operation which prohibits the use of sharply tuned passive 
devices. In this system the required isolation is obtained with a 
Hewlett-Packard Model C0043 3S03A microwave switch following the 
TWT power amplifier. This unit uses two PIN diodes in a common anode 
configuration and provides a 1 . 5 db maximum insertion loss and a 60 db 
minimum isolation across the 4-8 GHz band. 

The switch drive circuit used with the noise blocking switch is 
shown in Figure 4,4. Note that the circuit is basically the same as that 
of Figure 4.3, however the common anode configuration of this switch 
requires the use of PNP transistors in the output stage. 

4 . 2 . 1 . 7 Duplexer 

The duplexing or coupling of both transmitter and receiver to a 
common antenna is accomplished with a 3-part ferrite circulator. The 
circulator used is a Melabs Model H-317-264 which has a maximum 
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insertion loss of 0.4 db and a minimum isolation of 18 db across the 
4-8 GHz band. 

The relatively short minimum range (40 feet) of the system 
requires that the IF amplifier in the receiver not saturate during trans- 
mission since the recovery time is well beyond the minimum range. 

This obviously requires a much higher isolation during transmission 
than that provided by the ferrite circulator. As described previously, 
the high noise output of the CW TWT requires a high isolation during 
reception. The requirements for greater isolation than can be obtained 
with broad-band passive devices are met by switching both the trans- 
mitter and receiver lines. 

The complete modulation and duplexing arrangement including 
isolation switches are shown in Figure 4.5. The timing sequence is 
likewise shown in Figure 4.5. During transmission the combined isola- 
tion of the ferrite circulator and the receiver blanking switch is 94 db 
minimum. During reception the combined isolation of the ferrite circula- 
tor and the noise blocking switch is 78 db minimum. The need for 
relatively fast rise and fall times of the microwave switches and precise 
adjustment of the driving waveforms can be seen from the timing sequence, 

4 . 2 . 1 . 8 Antenna 

The antennas used are relatively standard 3-foot or 6-foot 
parabolic reflectors with broad-band (3,95 - 8.2 GHz), dual polariza- 
tion ridged waveguide feeds. The units are Scientific Atlanta Models 
22-6/A and 22-3/A reflectors with Model 28-3.9/6 and 28-3.9/3 feeds. 
E-plane patterns at the band extremes are shown for both antennas in 
Figures 4,6 and 4.7. 

The antenna is positioned for scattering measurements and 
scanned for imaging by means of a standard antenna positioner, Scientific- 
Atlanta Model 5403-1 . The positioner control. Scientific -Atlanta Model 
4112, and indicator assembly, Scientific-Atlanta Model 4 4 2 2 — (3 ) (3), 
permit a readout precision of +0.03° in both azimuth and elevation. 
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Modulation Noise Blocking 

Switch Switch 



Switch 



Antenna 



Figure 4.5. Duplexing arrangement and timing sequence. 




Angle 

(a) Frequency 4 GHz. 



Figure 4.6. E-plane patterns for 3-foot reflector. 
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Angle 

(a) Frequency 4 GHz, 



Figure 4.7. E-plane patterns for 6-foot reflector. 
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The positioner was modified for inverted mounting to conform to the 
geometry of the measurements performed. In this configuration the 
elevation axis permits viewing the ground from the nadir to grazing. 

The azimuth axis permits positioning and is also the scan axis for 
production of B-scan images. 

4.2.2 Receiver Section 

The receiver section consists of a low noise figure TWT pre- 
amplifier, microwave receiver blanking switch, bandpass filter, mixer- 
preamplifier, IF postamplifier-detector, and signal attenuators . 

4 . 2 . 2 . 1 Low-Noise Preamplifier 

The received signal is preamplified with a Watkins -Johnson 
low-noise travelling wave tube amplifier, Model WJ-423. This unit 
has a maximum noise figure of 7.3 db across the 4-8 GHz operating 
range. The minimum small signal gain is 31.8 db and the saturation 
power output is +13 dbm. 

4 . 2 . 2 . 2 Microwave Receiver Blanking Switch 

The minimum range capability of the system requires that the 
receiver be prevented from saturating during transmission. This requires 
a transmitter-to-receiver isolation considerably in excess of that provided 
by the duplexing circulator. Additional isolation is provided by a micro- 
wave PIN diode switch in the receiver line. The switch used is a series 
arrangement of two Hewlett-Packard Model 3570 microwave switches and 
drivers identical to those used for the modulation switch. 

4. 2. 2. 3 Bandpass Filter 

A bandpass filter is used prior to the mixer-preamplifier to 
eliminate video transients introduced onto the RF line by the microwave 
switching circuits. The filter is a Hewlett-Packard Model 8435A 4-8 GHz 
bandpass filter. This unit has a maximum passband insertion loss of 2 db 
and a rejection of over 45 db ten percent away from the band edges. 


81 


4.2 .2 .4 Mixer-Preamplifier 

The bandwidth of the IF amplifier must be great enough to 
accommodate the maximum bandwidth of the individual pulses. The 
minimum pulse length, , of 20 nanoseconds requires. approximately 
1 ,2/y or a 60 MHz bandwidth. In order to allow for possible future 
reduction of the pulse width an IF bandwidth of 100 MHz was used for 
the system. 

The incoming signal is translated to a center frequency of 250 
MHz. This mixing operation is accomplished in a coaxial integrated 
mixer-IF amplifier, RHG Electronics Laboratory Model MP 4-8/437. 

This unit has a signal bandwidth of 4-8 GHz, an IF bandwidth of 100 MHz, 
an IF noise figure of 4 db and a combined mixer-preamplifier power gain 
(RF to IF) of 23 db. 

4 . 2 . 2 . 5 IF Postamplifier- Detector 

The final stage of amplification is obtained with either a linear 
or logarithmic IF postamplifier. For most scattering coefficient measure- 
ments the greater precision linear amplifier is used. The mean of the 
fading return is positioned near the middle of the usable dynamic range 
by means of signal attenuators in both RF and IF lines. For most imaging 
applications the wide variety of returns as well as the range variation 
necessitates a much wider dynamic range. For this reason, a logarithmic 
final IF amplifier is provided for this application. 

The linear section is an RHG Electronics Laboratory Model EBT 
106 MGC IF postamplifier. This unit has a 250 MHz center frequency 
and a 100 MHz bandwidth . The amplifier gain is 55 db with a noise 
figure of 5 db. Both IF and video outputs are available with output 
capabilities of 1 milliwatt into 50 ohms and 2 volts into 90 ohms, 
respectively. 

The logarithmic section is an RHG Electronics Laboratory Model 
LT16040RFI IT postamplifier. This unit has a center frequency of 160 MHz 
and a bandwidth of 40 MHz, Amplifier linear gain is 62 db with a noise 
figure of 8 db. Video output is 2 volts into a 90 ohm load. When using the 
logarithmic amplifier, the 160 MHz IF frequency is obtained with a Sage 
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Laboratories Model 2543ND mixer. This unit has a maximum noise figure 
of 8,5 db and a maximum conversion loss of 8.5 db, 

4 . 2 . 2 . 6 Signal Attenuators 

The dynamic range required in making backscatter measurements 
over a range of incidence angles is far outside the capability of a linear 
system. Reference to the measured results presented in Chapter 2 shows 
that the scattering coefficient may easily vary from near +20 db at vertical 
incidence to less than -30 db at 60° incidence angle. Range variations in 
the measurements will add at least 10 db to this figure requiring a total 
dynamic range of approximately 60 db to measure differences in the mean 
value of the return. The dynamic range must accommodate not only the 
variation in c° but also the fading of the individual returns about this 
value. The 5% to 90% range assuming a Rayleigh distribution for the 
return is 18 db (+8 db to -10 db' . The total dynamic range required is 
in the order of 68 db. 

The system cannot handle this range even if the less accurate 
logarithmic IF postamplifier is used. In order to accommodate this range 
it is necessary to effectively step the receiver gain for different measure- 
ments . This is accomplished by placing variable attenuators in. the 
system prior to each amplification stage. Sufficient attenuation is placed 
in the system to prevent saturation of any amplifier and to place the mean 
signal level in the center of the available dynamic range. 

4,2.3 Local Oscillator Section 

In a swept frequency system such as this the local oscillator 

frequency must accurately track that of the transmitter to keep the 

heterodyned return signal within the IF passband. This is accomplished 

by developing the local oscillator signal from the same sources as the 

transmitter. A sample of the transmitted signal, f , is mixed with a 250 

MHz reference signal to produce sidebands at f + 250 MHz. This mixer 

output is passed through a narrow band tunable filter which tracts the 

upper sideband, f + 250 MHz. This offset sample of the transmitter 

c 

signal is then mixed with the return to produce the 250 MHz IF frequency 
used for subsequent amplification. 
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The local oscillator section consists of a PIN diode leveling 
loop, first mixer, reference signal source, and tunable filter-driver. 

4 . 2 . 3 . 1 PIN Diode Leveling Loop 

The sample of the transmitter signal is leveled using an 
independent PIN diode leveling loop identical to that in the transmitter. 
The reference signal used for leveling is taken off with a -10 db direc- 
tional coupler after the signal has passed through the tunable filter, 

4 . 2 . 3 . 2 First Mixer 

The sample of the transmitted signal is mixed with a 250 MHz 
reference signal using an Aertech Model 8000 mixer/modulator. This 
unit uses doubly balanced Schottky diodes and produces approximately 
1 milliwatt of power in the upper sideband across a signal input range 
of 4-8 GHz. 

4. 2. 3. 3 Reference Signal Source 

The 250 MHz reference signal is obtained from a Frequency 
Sources Inc. Model FS-7R cavity controlled mechanically tunable 
oscillator. Power output is approximately 200 milliwatts at 250 MHz 
with a 5% tuning range and a + 5% frequency stability over the operating 
temperature range of -30° to +70°C. 

4 . 2 . 3 . 4 Tunable Filter-Driver 

The upper sideband output of the first mixer is isolated by 
means of a four-stage YIG tuned filter, Watkins-Johnson Model ‘WJ- 
617-8. This unit has an integrated voltage-to-current driver, WJ-922-4, 
and tunes the frequency range of 4,25-8.25 GHz, The 3 db bandwidth 
is approximately 50 MHz with a maximum passband insertion loss of 5 db, 
an off-resonance isolation of 70 db, and a selectivity of 24 db per 
octave. The tuning voltage is 0-10 volts with a maximum deviation 
from linear of + 8 MHz. 
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4 . 2.4 Digital Control Section 


The frequency of the transmitted signal and the signal process 
sequencing are controlled digitally, A block diagram of the digital con- 
trol system is shown in Figure 4,8. 

A pulse generator, Hewlett-Packard Model 220A, serves as the 
master clock for the entire system. The generator is synchronized with 
the sweep frequency oscillator at e period of 1/1000 of the sweep 
period; thus, 1000 pulses are available for modulation of the transmitted 
signal. If the sweep is over the full 4-8 GHz range, the pulses are 
spaced 4 MHz apart. The frequency, f, of the transmitted signal may 
be selected by triggering the modulating PIN diode switch only after n 
pulses have been counted from the start of the sweep frequency period, 
Thus , 


f = 4 + — G Hi 

/ oo o 


(4.1) 


The effective bandwidth of the signal may likewise be deter- 
mined by counting the number, b, of consecutive pulses transmitted 
during a single sweep. 


MHZ 


(4.2) 


where BW 


BW 


eff 


bandwidth of the individual pulse 
effective bandwidth obtained by averaging of pulse 
sequence . 


Additional counters are used to automate data gathering by scanning a set 
effective bandwidth across the entire range of 4-8 GHz. 

The digital control modulation signal is fed to a pulse forming net- 
work that generates the driving pulses for all microwave switching. As 
considerable cable length may be present between the digital control 
unit and the microwave switch drivers, it is necessary to compensate for 
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pulse distortion effects which may occur. This is done by regenerating 
a fast-rise and fall time variable width pulse from the sharpened negative- 
going edge of the input signal. The circuit of this pulse forming network 
is shown in Figure 4.9. 

In addition, the digital control unit provides the sequencing 
signals for sampling of the return signal, the reading of the integrated 
average by a digital voltmeter and the drive signals and conversion 
circuitry required to record the digital voltmeter output on punched 
paper tape . 

4.2.5 Scattering Coefficient Data Processor 

In order to determine the differential backscattering coefficient, 

0° , it is necessary to measure the average backscattered power from the 
target area of interest. This is done by sampling the return signal at a 
time delay corresponding to the range of the target element. Signal 
averaging must be performed both to obtain the desired effective band- 
width and to eliminate effects of fading due to spatial variability of the 
individual scatterers within a resolution cell such as might be caused by 
wind motion. The averaging for increased effective bandwidth is performed 
on successive pulses within a single sweep while averaging of spatial 
fading effects is performed on a sweep-to-sweep basis. 

The functional block diagram of the scattering coefficient signal 
processor is shown in Figure 4.10. This unit is composed of a sampling 
oscilloscope, integrator, digital voltmeter, paper tape punch, and inter- 
face equipment. The signal process timing sequence is shown in Figure 
4.11. 

4. 2. 5.1 Sampling Oscilloscope 

A Tektronix Type 3S76 plug-in is used as the sample-hold circuit 

for the system. This unit has a sample gate width of approximately 0.6 

5 

nanoseconds and can operate at a repetition rate near 10 samples per 
second. The extremely short gate width insures that the measured 
amplitude is representative of a single independent sample of the return 
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V „ = -3 volts 
on 

V Qff = -2 volts 


All resistances in ohms 

Transistors: NPN 2N4420, PNP 2N4423 

Gates: Texas Instruments SN7410N 


Figure 4.9. Circuit diagram of pulse forming network. 
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Figure 4.10. Block diagram of scattering coefficient signal processor. 
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(a) Sampling sequence . 

Figure 4.11. 


Integrator React Command 

(b) Integration and read sequence. 
Signal process timing sequence. 
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distribution even when the individual pulse bandwidth is increased to 
occupy the full IF bandwidth of 100 MHz. The memory circuit of the 
pluc;-in retains the sample amplitude until the next sample period, a 
duration of approximately 10 microseconds at the maximum sample rate. 
This vertical signal sampled output is available at the front panel of 
the plug-in. 

The sampling plug-in is operated in conjunction with a Type 
3577 sampling sweep plug-in with the two units mounted in a standard 
main frame Type RM561 . Operated with an external trigger and in the 
manual sweep mode it is possible to manually adjust the time delay 
between the trigger and decayed sample trigger. This variable delay 
permits positioning of the sample gate at any desired position on the 
return signal . 

4. 2. 5. 2 Integrator 

Averaging of the return signal is accomplished by linear integra- 
tion. With a fixed repetition rate the sample width is constant, hence 
the mean value of the return may be determined by knowledge of the 
number of samples integrated. 

This application calls for an extremely stable integrator with 
provision for zeroing out any input offset voltage. Reset is accomplished 
by ci relay closure . The circuit diagram of the integrator is shown in 
Figure 4.12. 

4. 2. 5. 3 Digital Voltmeter 

A Hewlett-Packard Model 3440A digital voltmeter with plug-in 
Model 3443A is used to read the integrator output. This unit has a BCD 
recorder output that is used to convert the data to a paper tape output. 

4. 2. 5. 4 Digital Voltmeter-Tape Punch Interface 

The interface between the digital voltmeter and the tape punch 
provides code translation, level shifting and timing. The BCD output of 
the digital voltmeter is in a 1-2 -2 -4 code format. In order to have the 
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Operational amplifier is Fairchild 
Linear Integrated Circuit yA741. 


Figure 4.12 


SI 


paper tape output compatible with existing computer input devices, it 
is desirable that this code be translated into a standard ASC II code. 

The BCD output of the digital voltmeter has logic voltage levels 
of -1 volt and -24 volts. The Texas Instruments TTL integrated circuits 
used for all system logic functions operate with levels of 0 volts and -5 
volts, while the tape punch requires input levels of 0 volts and -10 volts. 
Obviously, level shifting is necessary if these units are to be operated 
in conjunction. 

The maximum rate of the tape punch is 12 characters per second 
and the system sample and read rates must not exceed this limit. The 
output readings are punched as 4 characters and a multivibrator clock is 
used to sequence the 4 characters to the appropriate punch driver at the 
maximum rate. 

All interface circuitry is included in the digital control unit 
shown previously in Figure 4.8. The circuitry described above is part 
of the logic I and logic II blocks shown in this figure . 

4.2 .5,5 Paper Tape Punch 

The punch used is a Commercial Control Corp Model N255 with 
a University of Kansas Model TPD/100 tape punch driver. The maximum 
rate of the punch-driver combination is 12 characters per second. 

It should be pointed out that paper tape is used as an output 
medium because of funding limitations which necessitated use of existing 
equipment. 

4.2,6 Image Data Processor 

The basic design of the system concentrated on obtaining back- 
scatter power spectra as a function of transmitted bandwidth and inci- 
dence angle. However, the applications of these data are principally 
to imaging radar systems. For this reason, an imaging capability was 
incorporated into the system design to facilitate demonstration of the 
techniques proposed. The images produced are uncalibrated and thus 
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must be judged on strictly a qualitative basis; however, the results 
shown in Chapter 6 are at least as striking as the earlier comparison 
of illumination by monochromatic and panchromatic light. 

The system records a B-scan presentation of the return signal 
as a Tektronix 585 oscilloscope. The cathode ray tube of this oscillo- 
scope has a special P-16 phosphor with rapid decay time and broad 
dynamic range. The image is recorded photographically directly from 
the cathode ray tube using a sheet film and developing process selected 
for maximum gray scale rendition. 

In normal imaging operation the logarithmic IF postamplifier is 
used to increase the usable dynamic range of the return signal. The 
detected video output is fed to a Tektronix 545 B oscilloscope operating 
in the delayed sweep mode. The vertical signal output of the Tektronix 
54 5B is applied to the cathode of the Tektronix 585 cathode ray tube to 
intensity modulate the beam. The delayed sweep of the Tektronix 545B 
permits positioning of the return signal so that the intensity modulated 
raster shows only the desired target return interval. The component 
arrangement and timing sequence for imaging operation are shown in 
Figure 4,13, 

4.3. System Calibration 
4.3.1 Calibration Technique 

The system is calibrated by referring all measurements to those 
made on a target of accurately known radar cross-section. The calibra- 
tion target used for this system is a 12-inch diameter polished aluminum 
sphere. A spherical target has the distinct advantage that its cross- 
section is not aspect-senstive since it remains a sphere no matter from 
what aspect it is viewed. The normalized radar cross-section of a 
sphere with infinite conductivity is shown as a function of relative size 
and wavelength in Figure 4.14 (Rheinstein, 1968). 
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Figure 4.13 Component arrangement and timing sequence for 
imaging operation. 
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Figure 4.14 Backscatter cross~section of an infinitely conducting 
sphere (after Rheinstein, 1968) . 
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The ultimate measurement desired is the backscattering cross- 
section, O', for discrete targets or the differential backscattering 
cross-section, 0°, for area extensive targets (Moore, 1969). The 
actual measurement performed by the system is of the average returned 
power, which is related to either O' or 0° by the radar equation. 

When measuring most discrete targets it may be assumed that 
the antenna gain is uniform over the target area and equal to the maximum 
gain; thus the familiar form of the radar equation may be used to 
give 


G* M = 


A* 


(4.3) 


where (/^.(A) 


p rt w 

P t (A) 

R 2 

A 

G(A,0) 


backscattering cross-section of the target 

return power 

transmitted power 

range from antenna to target 

wavelength 

antenna gain 


A similar measurement made on the calibration sphere yields 


0 S 60 


( 7?) 3 [. Pps Ca } f P t CO J P* 

A z G 2 (a, 0 ) 


(4.4) 


If we assume that during the course of the two measurements the trans- 
mitted power, receiver gain and antenna remain the same, equations 4.3 
and 4.4 may be combined to give 
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The radar backscattering cross-section of a discrete target may 
thus be referenced to the known backscattering cross-section of the 
calibration sphere without depending upon determination of the absolute 
value of each parameter in the radar equation. 

The measurement of the differential backscattering cross-section, 
<7°, of an area extensive target is slightly more complex as the antenna 
gain may no longer be considered a function of wavelength alone. Use 
of the differential backscattering cross-section implies that the return 
from the area extensive target is contributed by a large number of scatter- 
ing elements where phases are independent. If this condition is met the 
average return power may be computed by superposition of the power from 
the individual elements. If we assume that the scattering centers are 
sufficiently dense within the illuminated area we may pass from the 
summation to an integral over the illuminated area, thus 


P h (k & ) 



(47r) 3 R + 


d A 


(4.6) 


where (A, 0) = 
P T (A) 

G {AM = 
A 

(7° (A,6) = 

R 

A 


average returned power over an ensemble of surfaces 
transmitted power 

gain function of the antenna (assumed circularly 
symmetric with the angle from bore slight) 
wavelength 

differential backscattering cross-section 
range from antenna to target element 
area illuminated 


Technically would be noted as an average since this is a requirement 
for the superposition of power to be applicable, however, this designa- 
tion is dropped in line with normal convention and the average is 
merely implied. 
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Assuming that targets are selected such that (7° may be considered 
constant over the illuminated area, we have 


P h (\ c 
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(4.7) 


The return from the calibration sphere may be rewritten as 
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Combining equations 4.7 and 4.8 we have 


7°(\ e) = 


f,^) P M] 

«Sff^~- 4A 


(4.9) 


where g(y') is the normalized antenna gain function having a maximum of 
unity. The integral accounts for the "aperture effect" which is the 
weighting of the signals from the individual elements by the antenna gain 
function. This effect has been treated in detail by Parkins (1966). The 
"aperture effect" integral is calculated separately and is a constant for 
each A and © combination. 

Thus, the differential backscattering cross-section may like- 
wise be referenced to the calibration sphere without depending upon an 
absolute measure of the radar equation parameters. 


4.3.2 Accuracy of Measurements 

The fact that the backscattering coefficient or cross-section is 
determined by a ratio of the average target return and the calibration 
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return rather than an absolute measure of the return power considerably 
simplifies the calibration procedure. For a single frequency measurement 
with both target return and calibration return near the same amplitude and 
measured near the same time, the accuracy is estimated to be within + 1 
db. For swept frequency measurements whore the return amplitude must 
also be adjusted by attenuator insertion the accuracy decreases to 
approximately + 2 db. 

The principal cause of error is due to the method of frequency 
determination and the variation of system gain with frequency. Frequency 
is determined by subdividing the total sweep range (4 GHz) into 1000 
pulse intervals and then counting the number of pulses from the start 
of the sweep. The pulse period may drift between 999 and 1001 counts 
giving a frequency accuracy of approximately + 4 MHz near the end of 
the sweep. The setting of the start and stop points for the sweep intro- 
duces an additional error giving an overall accuracy for frequency deter- 
mination of near + 10 MHz. The gain variations of the system with 
frequency are most easily seen in the swept frequency calibration return 
shown in Figure 4.15. 

Another cause of error is the limited dynamic range of the receiv- 
ing section. With the linear postamplifier the usable dynamic range is 
approximately 20 db which is barely adequate for the expected fading 
range from complex targets. For targets of substantially larger radar 
cross-section than the calibration sphere or targets at a significantly 
greater range than the calibration measurement the average signal return 
is adjusted to fall near the center of the dynamic range by insertion or 
removal of signal attenuators in either RF or IF lines. The tolerance 
upon the attenuator values gives a corresponding uncertainty in the 
stability of the system gain which is assumed constant. 

To minimize the effect of the error sources, the data runs for 
each target were repeated a number of times. 
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CHAPTER 5 


► 

SCATTERING EXPERIMENT 


The experiment described in this chapter was conducted with 
two purposes in mind: the first, to verify the reduction of return 
variance obtained by frequency averaging, and the second, to investi- 
gate the spectral response of a variety of targets after averaging out 
the effects of fading. 

The results presented here were obtained with the experimental 
system described in Chapter 4 with two notable exceptions. During the 
period of data gathering the sampling oscilloscope and the boom truck 
were inoperative and undergoing extensive repairs. The loss of the 
sampling scope necessitated manual recording of the data with the aid 
of a storage scope. The results were then punched on computer cards 
for subsequent normalization and averaging. While this procedure did 
not hamper the reliability of the measurements, it did seriously slow 
the data gathering and limit the number of measurements that could be 
made „ 

The loss of the boom truck destroyed the system mobility and 
severely limited the choice of targets. The measurements shown here 
were obtained by mounting the entire system atop a ten-story dormitory 
and measuring suitable targets within the range of the system. This 
fixed mode of operation eliminated any investigation of the effects of 
incidence angle. Illustration 5.1 shows the system located on the 
dormitory roof. 

5 . 1 Description of Targets 

The field of view from the dormitory roof is shown in Illustration 
5.2. The following targets were measured and the location of each is 
shown by corresponding number in Illustration 5.2. 
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Illustration 5.1 Roof top installation of rnGasuromsnt systsm. 
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1. Space Sciences Building — This target was selected 
as representative of a large complex cultural target 
of size corrr «t- aole to or exceeding that of the system 
resolution cell. The range to the near wall of the 
building is 385 meters and the angle of incidence to 
ground level is 85 degrees. At this range the azimuth 
resolution is approximately 20 meters, while the 
range resolution is adjustable with a minimum of 

3 meters . 

2 . Grass — The grass in this area is regularly mowed and 
has a uniform height of about 2-1/2 inches. This 
target is representative of a uniform area extensive 
target with a low height vegetation cover and is, 
fortunately, located at a range of 78 meters. The 
angle of incidence is 63 degrees and the system 
azimuth resolution is approximately 4 meters. 

3 . Car — This target is representative of a small cultural 
target of a size less than the system resolution cell. 

For this measurement the car was located approxi- 
mately in the center of the far parking lot and headed 
directly at the system as shown in Illustration 5.2, 

At this angle of incidence (65 degrees) the parking 
lot surface provides an essentially zero return back- 
ground. The range is 80 meters and the azimuth reso- 
lution approximately 4 meters 

4. Car — The target and location are identical to that of 
target 3 with the exception that the car orientation is 
at right angles to the system line-of-sight. 

5. Area of Parked Cars — This target consists of a suf- 
ficient number of cars to totally fill the resolution cell. 

All cars are oriented to point approximately at the system. 
The angle of incidence, range and azimuth resolution are 
respectively 69 degrees, 97 meters and 5 meters. 
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6. 


Trees — This area is representative of a natural vegeta- 
tion target with a considerable vegetation cover. 

Unfortunately, the angle of incidence is 79 degrees 
which is slightly larger than the maximum incidence 
angle of most imaging systems. The range to the area 
is 180 meters and the system azimuth resolution is 
approximately 9 meters. 

7, Cement — This area represents a smooth, near-specular 
reflective surface over the wavelength range of this 
system. The angle of incidence to the area measured 
is 48 degrees and the system sensitivity is barely 
adequate to obtain the measurement at this angle. 

The cement parking lot surfaces located at a greater 
range and higher incidence angle are totally submerged 
in the system noise. The range to the area measured is 
53 meters and the azimuth resolution is approximately 
3 meters . 

5.2 M easurement Results 

The data recording was the same for each target with the only 
difference that of pulse length or actual resolution bandwidth. Voltage 
return measurements were recorded at 1000 points across the 4-8 GHz 
frequency range. The data were then normalized to the return from the 
calibration sphere and averaged for the desired effective bandwidth. 

The degree of averaging and the presentation of results depends 
upon whether one is interested in the reduction of fading or in the spectral 
response characteristic that remains after fading effects are averaged out. 
These data are analyzed separately for investigation of both phenomena. 

5.2.1 Measurements of Return Distribution 

In an earlier acoustic modeling experiement. Rouse (1968) attempted 
to measure the probability density function of the scattered return with 
frequency fading. The results of his experiment were inconclusive as the 
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frequency range of the acoustic system was inadequate to obtain sufficient 
independent samples of the frequency fading record to adequately define 
the distribution. With the limited range of his system the best that Rouse 
could do was to show that the addition of frequency diversity did not 
disturb the Rayleigh character of the positional fading distribution. 

With the increased frequency range of the radar system used for 
this experiment it was possible to obtain sufficient independent samples 
to reasonably define the frequency fading distribution through frequency- 
of-occurrence histograms. The probability density functions of the return 
distributions from each of the targets an; f hown in Figures 5.1 and 5.2. 

No averaging was performed on these data, thus the system resolution and 
effective bandwidth are equal. It may be seen from Figures 5.1 and 5.2 
that despite the dissimilarity of the targets, the behavior of the return 
distribution is much the same with each curve a reasonable approximation 
to the Rayleigh power or exponential density function. 

Of more significance to our investigations of variance reduction 
by frequency averaging is the behavior of the mean and standard devia- 
tion of the distributions. Table 5. T shows calculated values of the mean, 
standard deviation and their ratio for each of the targets measured. These 
figures were calculated from the data prior to averaging, thus the resolu- 
tion and effective bandwidth are the same. The theory of Chapter 3 pre- 
dicts that for a system with equal illumination and effective bandwidths, 
the ratio, (J-) , should be slightly less than one if the targets are truly 
representable by the model of Chapter 3 . It may be seen from the figures 
of Table 5.1 that this is not the case for these data. 

The explanation for this discrepancy must lie in either the system 
or the model. One might first suspect variations of the system sensitivity 
with frequency as the factor contributing to this increased variance; 
however, repeated measurements of the calibration sphere indicate the 
system response remains essentially constant with time. A more reason- 
able explanation is that the model makes no allowance for variations of 
scattering cross-section with large changes in frequency. Multiband 




(a) Target 5 - Area of cars 


(b) Target 6 - Trees 
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| (c) Target 7 - Cement 


Figure 5.2 Probability density functions of frequency fading distribu 
tions for targets 5-7. 




TABLE 5 1 

EXPERIMENTAL VALUES OF RETURN MEAN AND VARIANCE 


Target 

M 

<7 


1 

2.166 

3.206 

2.19 

2 

1,082 

2.104 

3.78 

3 

0.537 

0.630 

1.37 

4 

3.944 

5.314 

1.82 

5 

2.471 

6.290 

6.48 

6 

3,536 

5.062 

2.06 

7 

0.953 

1.624 

2.91 


data by previous experimenters Show, however, that such variations do 
exist, and the measurement range of this system approaches that of the 
discrete frequency measurements made in separate bands (Cosgriff et al., 
1960; Lundien, 1966). This argument is treated in greater detail in 
Section 5.3. 

The most graphic illustration of the effects of frequency averaging 
is examination of the full spectral response curve plotted for differing 
effective bandwidths. Figures 5,3 and 5,4 show the spectral response 
curves of targets 2 (grass) and 3 (car) plotted for effective bandwidths of 
100 MHz and 500 MHz. It may be seen from these plots that the slightly 
greater range extent of target 2 gives a higher fading rate and more rapid 
decorrelation with frequency shift, 

5.2,2 Measurements of Spectral Response 

Any measurement of scattering coefficient must include sufficient 
averaging to eliminate fading effects . In most single frequency systems. 
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Figure 5.3 Spectral response curves for two different effective bandwidths. 
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the fading is caused by positional changes due to relative motion between 
; radar and target. Averaging is performed in time with the assumption that 

during the averaging interval the surface being measured remains statisti- 
cally homogeneous . 

The system used here is stationary as are the targets , thus the 
only relative motion is that due to wind movement. For the cultural 
targets this motion is negligible at the wavelengths used for the measure- 
ments. The required averaging for these measurements must then be per- 
'v formed in the frequency domain . 

The pulse widths used for these measurements are all in the 
f range of 20 nsec, to 40 nsec, corresponding roughly to a resolution 

I bandwidth range of 50 MHz to 25 MHz. As shown in Chapter 3 # the 

reduction in variance by frequency averaging is approximately linear 
with multiples of the resolution bandwidth. Ten resolution bandwidths 
or 500 MHz was selected as the minimum effective bandwidth required 
to provide sufficient averaging of the data to remove fading effects and 
to permit viewing of the underlying spectral response. 

Figures 5.5 through 5.11 show the spectral response character- 
istics of the seven targets described. If there is a single dominant 
characteristic of these data it is that, in general, the scattering coefficient 
increases with increasing frequency. 

1 

5.3 Discussion of Results 

In this section an attempt is made to relate the experimental 
results with the theoretical predictions of Chapter 3. This, of course, 
applies only to the investigation of panchromatic illumination or the 
reduction of return variance with frequency averaging. No theory is 
put forth to explain the spectral response characteristic of such a diverse 
array of targets and indeed, it is doubtful that this is possible. The 
investigation of spectral response was rather an attempt to assay the 
potential value of polypanchromatic illumination, or the application of 
color techniques, to the microwave spectrum. 
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Figure 5 . 5 Spectral response of target 1 
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Figure 5.7 Spectral response of target 3. 



Target: Car (Broadside) 
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Figure S.10 Spectral response of target 6 
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Figure 5.11 Spectral 
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of target 7. 







5.3.1 Variance Reduction with Panchromatic Illumination 


Frequency averaging will reduce the variance of the return signal, 
but Hie question remains: how much bandwidth is required for a given 
reduction of return variance and how well do the results and model of 
Chapter 3 predict the observed behavior? 

The discrete scatterer model proposed in Chapter 3 was shown 
to give rise to a Rayleigh distributed return when monochromatically 
illuminated. The probability density functions shown in Figures 5.1 and 
5.2 seem to indicate that the Rayleigh distribution is at least a reason- 
able approximation for the targets measured. However, the mean and 
variance figures shown in Table 5.1 seem to contradict this appearance 
for a Rayleigh power distribution has equal mean and variance. In 
addition, it was pointed out in Chapter 3 that a system with a resolving 
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capability should actually have a ratio, (—) , of slightly less than one. 

The answer to this contradiction lies in the frequency variation 
of the scattering cross-section noted in the spectral response curves. 
Figures 5.5 through 5.11. The calculations of mean and variance shown 
in Table 5 . 1 treat the data as a function of a single variable which is 
clearly not the case. Rather, the data represents at least a two-variable 
distribution where the variables are large and small scale variations of 
frequency. The distribution considered in Chapter 3 arises from relatively 
small changes in frequency which affect the relative phase between 
scattering centers in the resolution cell. However, no allowance was 
made for a change of scattering coefficient that may accompany frequency 
changes large enough to substantially alter the reradiation pattern and 
gain of the individual scattering elements. 

To investigate the reduction of fading variance alone we must 
modify our calculation of variance to account for larger scale frequency 
variations . It would appear this might be done by fitting the data with 
a regression line and considering the variance about this regression line 
to represent the fading variation predicted by the theory of Chapter 3 . 


120 


There still remains the choice of the type of regression analysis which 
should be used. From the spectral response curves it appears that at 
least two scales of variation are occurring above that of the fading. On 
a scale of a few hundred megahertz or less there appear to be appreciable 
resonance effects even with an effective bandwidth of 500 MHz. On an 
even broader frequency scale of at least a few gigahertz there appears a 
variation that might be attributed to the effective size of the scattering 
elements . 

While a simple linear regression may suffice for those curves 
with no noticeable resonance effects, others will require a polynomial 
fit of the data. The selection of the degree of the polynomial is some- 
what arbitrary and will directly affect the magnitude of the variance 
attributed to fading eifects . 

As an alternative to performing the regression analysis the 
variance was recalculated using an effective bandwidth curve of 1 GHz 
as an approximate regression line. The resulting mean, variance and 
their ratio are shown in Table 5.2. 

The figures shown in Table 5.2 agree much better with those 
predicted by the analysis of Chapter 3. As expected, the greatest 
improvement is realized for those targets with the steepest slope in the 
spectral response characteristic. The ratios, are still somewhat 

larger than predicted; however, it is felt that this is due to differences 
between the resonance scale variation and the still larger scale approxi- 
mately linear variation. The use of an effective bandwidth large enough 
to remove all resonance range effects gives approximately the same 
variance as the figures shown in Table 5.1. The use of a smaller 
effective bandwidth would undoubtedly yield better correlation for the 
figures of Table 5.2, however to go below 500 MHz would yield results 
of questionable validity as the approximate regression curve would 
actually be following to an extent the fading variations . 

The improvement in return variance with increasing effective 
bandwidth is illustrated by the figures of Table 5. 3. The figures shown 
are the predicted and measured Improvement in (£)? between an effective 
bandwidth equal to the resolution bandwidth and the effective bandwidths 
shown. 



TABLE 5 . 2 

EXPERIMENTAL MEAN AND VARIANCE OF FADING DISTRIBUTION 
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Target 


or 

(Zy 

V 

1 

2.166 

3.028 

1.96 

2 

1.082 

1.294 

1.43 

3 

0.537 

0.617 

1.32 

4 

3.944 

5.057 

1.64 

5 

2.471 

1.992 

0.65 

6 

3.536 

4.567 

1.67 

7 

0.953 

0.981 

1.06 


TABLE 5 . 3 

MEASURED AND PREDICTED IMPROVEMENT IN VARIANCE 
AS A FUNCTION OF ILLUMINATION BANDWIDTH 



Eff. BW i 

= 100 MHz 

Eff. BW = 

300 MHz 

Eff. BW 

= 500 MHz 
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/(“)? 
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/(-) 2 
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Target 

Predicted Measured 

Predicted Measured 

Predicted Measured 

1 

2.0 

1.5 

5.3 

3.9 

8.8 

6.1 

2 

2.0 

2.4 

5.3 

6.4 

8.8 

9.2 

3 

2.0 

2.3 

5.3 

5.5 

8.8 

6 . 9 

4 

2.0 

to 

• 

o 

5.5 

5.0 

8 . 8 

6.5 

5 

2.0 

2.4 

5.3 

4.9 

8.8 

7.9 

6 

2.0 

1.9 

5.3 

3.9 

8.8 

5.3 

7 

2.0 

2.0 

5.3 

4.6 

8.8 

6.3 
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While the agreement of the predicted and measured values is not 
precise it is reasonably close. This is to an extent to be expected since 
the analysis of Chapter 3 shows almost any averaging to follow an approxi- 
mately linear decrease in 

The results of the experiment conclusively show the benefit of 
panchromatic illumination. The agreement with the theory developed in 
Chapter 3 is not so apparent. This is felt to be a result of the additional 
frequency sensitivity ranges encountered rather than a deficiency of the 
model used to predict the fading variations . A truer test of both the value 
of panchromatic illumination and the applicability of the model would be 
provided with a system capable of absolute bandwidth variations over the 
range used and with the additional capability of positional averaging through 
relative movement. Such a system would be capable, at least for area 
extensive targets, of measuring an ensemble of homogeneous elements, 
thus reducing the resonance effects encountered here. 

5,3.2 Spectral Response Curves 

The investigation of the spectral response of radar target cross- 
section was undertaken in part to assay the potential signature content 
of such data. If recognition is to be attempted on the basis of these data 
the response must exhibit sufficient variation across the frequency range 
to discriminate between various targets. Examination of the spectral 
response curves of Figures 5,3 through 5.11 shows there to be distinct 
differences in the response characteristics of the different targets. 

Let us for the moment ignore the differences in mean level 
between the curves and concentrate on the shape of the response 
characteristic . It appears that the variations shown occur on three 
different frequency scales. First is the very fine scale fading variation 
discussed in the previous section. This is seen to still be quite prevalent 
in the curves for 100 MHz effective bandwidth. In the curves for 500 MHz 
effective bandwidth we may assume that the fine scale fading variations 
have ess mtially been removed by the averaging. We note, however. 
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that the variations remaining still have two distinct frequency scales. 
Apparently there are resonance effects that give pronounced variations 
on a scale of a few hundred megahertz. Superimposed on this resonance 
characteristic is still another variation on a scale of at least a few 
gigahertz. This latter trend was in the previous section attributed to 
the frequency sensitivity of the individual scattering elements. 

To attempt recognition on a frequency scale that still contains 
appreciable fading variation appears to be a hopeless task. Either of 
the two broader scale variations might be used in an attempt to define 
spectral signatures. However, the resonance region variation, on a 
scale of a few hundred megahertz , would in all likelihood be overly sensi- 
tive to slight changes in target configuration or orientation. In addition, 
calibration of this range will be much more difficult than that of the 
broader scale. 

It appears that recognition of spectral response characteristics 
would be both easier and more meaningful if the variations considered are 
on the scale of at least one gigahertz. To further investigate this scale 
of variation the spectral response curves of all targets were replotted 
with an effective bandwidth of 1 GHz. These curves are shown in Figures 
5.12 and 5.13. Differences between the targets measured are still 
readily apparent on this scale . 

Before proceeding further with possible recognition schemes for 
the spectral response characteristic let us review the most universal 
means of spectral response recognition, namely the perception of color 
in the visual region. The great variety of color as perceived by the human 
eye arises from spectral response variations across a region less than 
one octave wide . The receptor of this information is sensitive to a 
wavelength range of only 0.4 to 0.7 microns. 

The instrument most often used to detect visual range response varia- 
tions is the camera with color film. Colors are produced, in the case of 
film at least, by separating the response of the region into three response 
ranges which are appropriately dyed and weighted to produce a result 
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| Figure 5.12 Spectral response with 1 GHz effective bandwidth. 
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equivalent to that perceived by the eye. Figure 5.14 illustrates the 
inherent color sensitivities of a three -layer film and the effective color 
separation of the emulsion layers (Sorem, 1967). 

The type of response variation that is sensed by either film or 
eye is illustrated by the spectral reflectance curves of Figures 5,15 and 

5.16 (Gates, 1970; Holmes, 1970). At first glance it may appear that 
these curves have significantly more variation than the curves shown for 
the microwave region of the spectrum . However, it may be shown that 
this is principally due to differences in the conventional means of present- 
ing spectral reflectance and radar cross-section data. While reflectance 
is normally shown on a scale of 0 to 100 per cent, radar cross-section 
data is normally presented on a compressed logarithmic scale. The effect 
of this scale difference may be illustrated by replotting the radar cross- 
section data of Figures 5.12 and 5.13 in terms of normalized ratios similar 
to those shown for the visual region. These curves are shown in Figures 

5.17 and 5.18. The normalization of each curve was with respect to its 
own mean value across the range; thus, the absolute values between curves 
do not correspond and the curves should be compared in functional form 
only similar to those of Holmes (1970), Figure 5.16. 

From the above it may be seen that the potential variation in the 
microwave region is easily as great as that currently being exploited in 
the ultraviolet, visual and infrared ranges. This Is apparent from an 
investigation of only a single octave of the microwave spectrum. One may 
only guess what the potential may be if greater use is made of the micro- 
wave region with a range of perhaps ten octaves. The techniques of 
recognition applied in the visual region may easily be adapted for use 
in the microwave region as will be demonstrated in Chapter 6. 
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Figure 5,14 Color film dye layer sensitivity ranges 
(from Sorem, 1967). 
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(a) Spectral reflectance of seven species of plant 
representing diverse groups . 


Acariutmim (bark) 
Outran •«» (bark) 

Pika Fur (lowar) 

Stallar Jay Faathan (back) 
Sand Dollar 
Human Hand 
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(b) Spectral reflectance of bark from two species of 
trees and parts of several animals. 


Figure 5.15 Spectral response variations in the visual region 
(from Gates, 1970) . 
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RELATIVE SPECTRAL RADIANCE 



Figure 5.16 Relative spectral radiance spectra for agricultural 
scenes in the range of 0 . 4-1 . 05u , August 30, 1966, 
11:50-11:56 a.m. Curves should be compared in 
functional form only; gain settings were changed 
between spectra (from Holmes, 1970). 
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Figure 5.17 Normalized spectral response curves. 
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CHAPTER 6 


IMAGING EXPERIMENT 


The experimental results of Chapter 5 prove the original con- 
tention that frequency averaging may be used to reduce return variance 
and that spectral variations sufficient for recognition purposes are 
present in the microwave region. These results are not confined to 
any particular type of radar system and may find numerous applications. 
The application of these techniques to imaging radar systems was 
selected for demonstration as having great utility and giving a vivid 
and graphic presentation of the results. 

In recent years the use of imaging radar systems for military 
reconnaissance and geoscience investigations has expanded rapidly 
(MacDonald, 1969; Wing, 1970). The imaging system provides an 
excellent example of the image improvement afforded by frequency 
averaging. A fully focused monochromatic synthetic aperture system 
obtains a single sample of the positional fading distribution. This 
gives rise to a large return variance which is most striking in the case 
of the "salt and pepper" effect noted in images taken by fully focused 
synthetic aperture systems. This effect seriously hampers the utility 
of such systems for applications such as change detection or the de- 
lineation of homogeneous area. The addition of panchromatic illumination, 
or excess range bandwidth, yields additional independent samples of the 
return distribution and thus reduces return variance and the scintillation 
noted in imagery of even homogeneous area extensive targets . Panchro- 
matic illumination provides the radar imaging system designer another 
degree of freedom in his balancing of trade-offs between such parameters 
as swath width, power, resolution, image, quality, etc. 

The interpreter of radar imagery is interested in the 
recognition or identification of objects portrayed on the image. Much 
has been written regarding the use of a possible multispectral signa- 
ture for the purposes of this identification. In addition , a number of 
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current programs are investigating partitioning of the spectrum in the 
infrared, visible and ultraviolet regions as a means of accomplishing 
such recognition. The experiment conducted here was designed to show 
that an equivalent amount of signature information resides in the micro- 
wave region and that identification techniques such as the use of false 
color imagery are applicable here as well. 


6.1 Measurement Technique 

The method of obtaining images was briefly described in Chapter 4 , 
The return signal is detected, amplified with a broad-bandwidth oscillo- 
scope and applied to the cathode of another broad-bandwidth oscilloscope . 
The imaging oscilloscope is modulated with a horizontal deflection voltage 
proportional to the angular position of the azimuth axis of the antenna 
positioner while the vertical axis is modulated with a ramp voltage propor- 
tional to the system range. The large difference in sweep times produces 
a raster display with many essentially vertical range lines for each 
horizontal azimuth line . Each return pulse intensity modulates a single 
vertical line . The result is a conventional B-scan display where an angular 
illuminated sector is presented in a rectangular display with azimuth as 
the abscissa and range the ordinate. 

The PRF of the system is set at 100 kHz while the azimuth scan 
rate is approximately ten degrees per second. This large difference in 
rate gives a separation of the vertical range lines far below the resolution 
capability of the CRT or the recording film. This results in an effective 
integration of successive pulses on the CRT-film combination. When 
panchromatic illumination is desired the frequency is swept at a rate such 
that the entire sweep range is covered in less time than that required to 
traverse an azimuth resolution cell. The PRF of the system was held 
constant during the production of all images shown here, thus the sample 
rate and effective transmitted power were the same for monochromatic 
and panchromatic illumination. 
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A few limitations of the system imaging capability should be 
kept in mind when viewing the data. As mentioned previously* the 
presentation is a conventional B-scan which distorts an illuminated 
angular sector into a rectangular format. The range bandwidth of the 
system in imaging operation is constrained by the amplifying oscillo- 
scope to approximately 25 MHz maximum. The combination of a rela- 
tively low system altitude and the use of a conical beam antenna pre- 
vents imaging of any appreciable range 

near grazing. At lower incidence angles the vertical beamwidth of the 
antenna confines the illuminated sector to a very narrow range dimen- 
sion . With the above limitations on imaging operation it should be 
realized that the images shown here are not presented as a calibrated 
data format, but rather are intended merely to demonstrate the potential 
application of panchromatic and polypanchromatic illumination to air- 
borne imaging systems. 

The system was operated with a conical beam antenna of 
approximately 3.2° beamwidth at the center frequency. The pulse width 
was 2 5 nsec., however, the full resolution capability of this pulse was 
not realized due to the bandwidth limitation of the amplifying oscillo- 
scope, The depression angle used was approximately 10° which gives 
sufficient range illumination to provide reasonably extensive image 
coverage. 

The area imaged is that shown in Illustration 5.2. The angular 
coverage is approximately 90° with the near range of the image commenc- 
ing approximately at the highway and the far range extending to the 
space sciences building. 

6.2 Measurement Results 


The results clearly indicate the improvement of image quality and 
definition obtained by the addition of frequency averaging or panchromatic 
illumination. Figure 6.1 shows successive monochromatic images of the 
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area made with no change in system parameters . Close examination of 
these images indicates differences even though there is no relative 
motion between system and target other than the slight amount provided 
by the wind. For the purposes of orientation the bright line labeled "A" 
on the image is the far bank of the highway with the overpass visible in 
the near range. The target labeled “B" is the Space Sciences Building 
measured in Chapter 5 . The differences shown here with such limited 
movement indicate the difficulty of performing change detection or 
recognition with a monochromatic system. 

For the purposes of comparison. Figure 6.2 shows a sequence 
of monochromatic images taken with slight differences in frequency. 

The changes due to the relatively minor motion caused by wind are 
indistinguishable from those produced by small changes in frequency. 

Figures 6.3 through 6.5 illustrate the improvement possible with 
panchromatic illumination. Each figure shows a monochromatic image 
and two panchromatic images with effective bandwidths of 0.5 GHz and 
1.0 GHz. The improvement in image quality is so obvious as to hardly 
require comment. The images of Figure 6.3 are centered about 4.5 GHz 
while those of Figures 6.4 and 6.5 are centered at 6.0 GHz and 7.5 GHz, 
respectively. 

The three 1 GHz bandwidth images may be considered the three 
components of a polypanchroma tic system. As such it should be possible 
to combine the information on a single image and perceive differences 
due to variations in the spectral response of the various targets. While 
an effort has been made to compensate for variations of system response 
between the three bands, the images should not be considered as having 
an absolute reference of better than about 3 db . 

Figure 6. 6 shows a color combination of the three panchromatic 
images shown. While the calibration of this image is crude at best, 
the variation of colors noted is an indication of the response recognition 
possible with polypanchroma tic illumination. 

Figure 6.7 is a sketch of the area imaged showing the relative 
target locations for the distorted B-scan images. 


136 



Figure 6.1 Successive monochromatic images made with no 
change in system parameters. 
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(c) f = 6.05 GHz. 


Figure 6.2 Monochromatic images taken at slightly 
different frequencies. 
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(a) Monochromatic f= 4.50 GHz. 



(b) Panchromatic f = 4 . 50 + 0 . 25 GHz . 



(c) Panchromatic f = 4 . 50 + 0 . 50 GHz . 


Figure b.3 Comparison of monochromatic and panchromatic 
images . 
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(a) Monochromatic f= 6.00 GHz. 



(b) Panchromatic f = 6 . 00 + 0 . 25 GHz . 



Figure 6.4 Comparison of monochromatic and panchromatic 
images . 
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Figure 



(a) Monochromatic f = 7.50 GHz. 



(b) Panchromatic f = 7. 50 + 0.25 GHz. 



(c) Panchromatic f = 7 . 50 + 0 . 50 GHz . 

6.5 Comparison of monochromatic and panchromatic 
images . 
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Figure 6.6 Polypanchromatic images. 
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(a) Undistorted sketch of target area. 



(b) Sketch of target area in B-scan geometry. 
Figure 6.7 Relative target locations on B-scan imagery. 
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6.3 Discussion of Results 


The results shown are largely self-explanatory* The images shown 
were obtained with far from optimum conditions. Namely , the resolution 
and dynamic range of the recording oscilloscope were marginal while the 
limited height of the system required operation over excessively high 
incidence angles. Despite these handicaps the improvement in image 
que lity obtained with panchromatic illumination is quite striking . This is 
particularly evident in the recognition of linear features such as those 
shown in the parking lot just beyond the overpass. Additional increases 
of effective bandwidth are observed to improve the definition or fill/ 
particularly of large complex targets . 

The curves of spectral response variation shown in Chapter 5 
illustrate the point that the variations encountered in the microwave 
spectrum are easily as large as those encountered in the visual spectrum. 
The potential use of this variation to produce color is at once immediate 
and obvious . The polypanchromatic color imagery shown demonstrates 
the ease with which color or any other recognition techniques used in the 
visual region may be applied to the microwave spectrum. 


CHAPTER 7 


CONCLUSIONS 


The goal of this work was two-fold: first, to investigate the 
effects of frequency averaging upon the variance of the return signal; and 
second, to investigate the continuous spectral response of targets in 
the microwave portion of the spectrum and to compare this response with 
that observed in the visual region. 

The frequency response of several targets was measured over a 
4 GHz frequency range and with a sampling bandwidth of 40 MHz . The 
fine scale frequency \ariations (the frequency fading statistics) were 
measured and compared with those predicted by a theory based on a 
discrete isotropic scatterer model. The decrease of return variance with 
increasing illumination bandwidth was predicted quite well; however, the 
absolute value of the measured variance was higher than predicted for 
all targets . This discrepancy was attributed to additional resonance 
variations with broad scale frequency changes not taken into account 
by the theoretical model. The use of an approximate regression analysis 
for the calculation of the frequency fading variance substantially 
improved the agreement of measured and predicted absolute values . 

The reduction in variance was shown to be essentially the same 
for a frequency averaging system and one operating at maximum resolution 
with post-detection integration. The advantage of the frequency averaging 
or panchromatic illumination is not in more efficient variance reduction 
but in possible decreased receiver and recording bandwidths. 

The utility of frequency averaging was demonstrated by producing 
B-scan imagery of an area with both monochromatic and panchromatic 
illumination. The improvement in target and pattern definition obtained 
with the reduced- variance panchromatic illumination is obvious . 


Examination of the spectral response characteristics led to the 
conclusions that at least three scales of variation with frequency were 
occurring: the fading variations predicted by the theory of a random 
scatterer model * a scale of a few hundred megahertz attributed to resonance 
phenomena, and a gigahertz scale variation attributed to changes of the 
reradiation pattern of the basic scattering elements. 

Comparison of the observed microwave variations with those 
measured in the visual range revealed that the ranges of variations are 
comparable. The potential of color as a means of spectral response 
recognition is obvious . This is merely one example of a visual spectrum 
technique applied in the microwave region. The range of microwave 
variations shown indicate that virtually any multiparameter spectral 
recognition technique is applicable to the microwave spectrum . 

The experiment performed resulted in defining a number of areas 
requiring further work. The most important of these are the following: 

First and foremost is the need for additional continuous spectral 
response measurements in the microwave region. The measurement 
program should be expanded to cover a frequency range of at least 0.4 
GHz to 35 GHz and must include a wide variety of target types. 

Spectral response measurements over an ensemble of homogeneous 
target elements are needed to determine the effect upon the resonances 
noted. This will require a system with the additional capability of move- 
ment during measurement. 

The application of panchromatic averaging to imaging radar systems 
will require further definition of the averaging means to take best advan- 
tage of the bandwidth compression to reduce system requirements and 
complexity. 
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APPENDIX A 


DERIVATION OF AN IDENTITY 

In equation (3.30) we encountered an expression of the form 

£ I ( A. 1) 

where Z^ and Z ^ are complex gaussian variables with zero means and 
uniform phase distributions. 

By the use of the correlation coefficient, p, we may separate Z 2 
into two orthogonal components. 

2 * = pi, + 2 . 

where 

£[».«.* 1 

E[I2,I*] 

Since the process is gaussian, Z^ and Z 2 are statistically independent* 

Thus we have 

- 1 z,i z li.l z a.*.* +/»/z,/*a*2. (a. 2 ) 


* 
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and 




£ [n l i*ii'i'] =* p’eLu.i*] + eLui^eCh.i 1 ] 

(A. 3) 

As the distribution has zero mean, the cross terms drop out leaving 

£ [/i,i*iiJ*] zp'ELiu’l+Ebi.i'Ubij*] (A 4) 

It was shown in Chapter 3 that the envelope of a complex 
gaussian process with zero mean gave rise to the Rayleigh distribution 
function . 



(A. 5) 


The fourth central moment of the Rayleigh distribution may be expressed as 

£ [ X 4 ] c ZE*[* 2 ] (A. 6) 

Recall also that Z 0 = i z -f2 t 


We may now rewrite equation (A. 4) as 






(A. 7) 
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* 


Recalling also that 

p £•£/*, / l J = E[i,i?J 

and that eCi.t *]=£&* it] we may write 

E[ II, 1*1211*] - SE* [l,Z*] + E [I hi*] e[HzI z ] 

- * f '[*. it*] 

and finally sine?; f f / 2i/ J J ^0 hava 

* £[2 l J + E*[i,l,*] 


(A. 8) 


(A. 9) 
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APPENDIX B 


EVALUATION OF AN INTEGRAL 


Equation (3.46) is given by 

4 1FR Z 

r i 2<r c 4 r ***** 




TTR 


f 


* '* 


['- jk] 


(3.46) 


Consider the integral 


wR 


i - £ J ^ [>■ 4-J 

o 


(B.l) 


Let 




WR 

2 f 

TTR J ** 


irR 


n*R* /' 


S!h z X 




(b. 2) 


Then 


ttR 

T 2_ f J_-coslX 

1 i ^ 7TR J 


2X 




(B .3) 


Integrate by parts letting 

= / * c os ZX 


dvr = 


oUl 3 2 X/ /» 2* 

t 


I - jL [ 

*• ) TR L 


Z* z 


cos 2 * - I 


TT - - 


2 * 


)T 0 ^ 

r * ./i 


2* 


Sth Zx 


ZX 


(B .4) 


For 1 2 we have 



irR 

^ / cof ix - I 

w*n z J i* 
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J*K 


(B . 5) 


I is now 
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(B.6) 


Let p = 2x and Q * 2irR, then 



£c$,f Q ~ t ] + 



alii + 


Q 



COSJU - I 

7c 


du 

(B . 7) 


Rational approximations are used for the sine and cosine integrals 
taken from National Bureau of Standards AMS 55 , Handbook of Mathematical 
Functions with Formulas, Graphs, and Mathematical Tables. We have 
from here 


<9 


/ SthM . IT - 

* fm COS0 - J (Q)sthQ 


(B.8) 


Q 

Ci (Q) s \r + l*Q+ * { (R)sthR -j(Q) *°SQ 


9) 
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where y = 0.57721567 and 


f(0) * 



I f ff 4 + 7 -24//6J <?**«?. 4*397* 7 _ 

<7 I* G 4 t 9. 06 f rfo <? 2 + 7. 1*7433 J ^ ^ 

/ I" <? 4 + 7. *4 7 47? <3? 2 + /.<r**0 72 ] / 

l <?4 + /i“. TtTTol <? 2 + /2. 723**4 J * "cp 8 


(B.10) 


(B.ll) 


We have finally for I 


X = JL ( £,(<?) + -~r [ ci(Q) ~ r - /h o +cosq -/ ] 

Q ( <? L J 

= i H' i ">**+iz* ,n «]~ B t i* s,n9 *o* coi9 1 


+ i J + ri 


(B.12) 


The program used tc evaluate I for a range of values of R and Q 
is included. 


152 


1 DIMENSION R (41) 

2 READ(5, 1) (R(I) ,1-1,41) 

3 1 FORMAT (10 (F. 2)) 

4 WRITE (6, 3) 

5 3 F0RMAT(1H1,8X,1HR,11X,1HQ,15X,3HSIQ,17X,3HCI0,18X,4HDARI) 

6 DO 10 K=l,41 

7 Q~R(K)*6.28318531 

8 A=((Q**4)+7.241163*(Q**2)+2.463936)/ 

9 1 ((Q**4)+9, 068580* (Q**2)+7. 157443) 

10 B-((Q**4)+7. 547478*(Q**2)+1 *564072)/ 

11 1((Q**4)+15.723606*(Q**2)+12. 723684) 

12 SIQ-l . 57079633- (A/Q)*COS(Q)-(B/Q**2)*SIN (Q) 

13 CIQ=(A/Q)*SIN(Q)-(B/Q**2)*COS(Q) 

14 DAIR=VARI (Q , A , B) 

1 5 WRITE (6 , 2)R(K ) , Q ,SIQ, CIQ ,DARI 

16 2 FORMAT(5X,F7.2,5X / F8.3,5X # E15.8 # 5X # E15.8 / 5X / E15.8) 

17 10 CONTINUE 

18 STOP 

19 END 


1 FUNCTION VARI(P,A,B) 

2 VARI=(4,O/P)*(l,57059633+A*((SIN(P))/(P**2)-(COS(P))/P)-B*((SIN(P) 

3 l)/(P**2)+(COS(P))/(P**3))+(l. 0/P)*(COS(P)-1 ,57721567-ALOG(P))) 

4 RETURN 

5 END 
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